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Abstract 

We study the boundary correlation functions in Liouville theory and in solvable statistical 
models of 2D quantum gravity. In Liouville theory we derive functional identities for all 
fundamental boundary structure constants, similar to the one obtained for the boundary two- 
point function by Fateev, Zamolodchikov and Zamolodchikov. All these functional identities 
can be written as difference equations with respect to one of the boundary parameters. 
Then we switch to the microscopic realization of 2D quantum gravity as a height model on 
a dynamically triangulated disc and consider the boundary correlation functions of electric, 
magnetic and twist operators. By cutting open the sum over surfaces along a domain wall, we 
derive difference equations identical to those obtained in Liouville theory. We conclude that 
there is a complete agreement between the predictions of Liouville theory and the discrete 
approach. 
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1 Introduction 



The two-dimensional quantum gravity (2D QG) allows two complementary descriptions based 
on Liouville field theory or on dynamical triangulations [2]. Prom the point of view of string 
theory, 2D QG is the theory of non-critical bosonic strings, which do not have a tachyon state 
and therefore admit an interpretation as statistical ensemble of embedded worldsheets. (For 
reviews see (HUH) In this paper we will consider the boundary problem in the two descriptions. 
We will present a systematic approach to evaluate the boundary correlation functions in Liou- 
ville and discrete 2D QG based on a set of difference equations with respect to the boundary 
parameters. Although derived by completely different methods, the difference equations for the 
two- and three-point functions obtained in both approaches become identical after rescaling with 
an appropriate factor. This factor does not depend on the Liouville boundary parameters and 
has a natural interpretation as the contribution of the matter fields. 

1.1 The continuum approach 

The first success of the continuum approach was to predict the scaling exponents of gravita- 
tionally dressed primary fields OEI • Subsequently, the so-called resonance correlation functions 
of vertex operators were found by doing the zero mode integration over the Liouville field and 
then using analytic continuation with respect of the number of the Liouville interaction inser- 
tions 0IB1- These results showed that the Liouville CFT, although 'noncompact', can be treated 
similarly to the standard 'compact' conformal field theories (§]• Indeed, in a series of impressive 
works |10|I11|I12|I13|I14| . the Liouville CFT has been practically solved. The exact expressions ob- 
tained for the bulk jlUUllj and boundary |12U13llT%| Liouville structure constants were obtained 
as solutions of a set of identities that follow from the OPE of vertex operators with degenerate 
bulk or boundary Liouville fields. 

In Liouville theory with boundary the situation is more subtle because of the presence, be- 
sides the cosmological constant fj,, of a second dimensional parameter, the boundary cosmological 
constant The boundary Liouville structure constants thus depend on a new dimensionless 
parameter, [ibI \JJi- The observables are meromorphic functions of this parameter having a 
branch point singularity at the point — fjfg where hb equals the boundary entropy density. The 
branch point can be resolved by introducing a uniformization parameter r 



In the correlation function of n boundary operators on the disc, one can define n indepen- 
dent boundary cosmological constants [abi, ^Bn associated with the segments between two 
neighboring operators. Fateev, A. Zamolodchikov and Al. Zamolodchikov ^2] found that the 
boundary two-point function, or boundary reflection amplitude, satisfies a simple functional 
equation, which can be given the form of a difference equation with respect to one of the bound- 
ary parameters, t\ or T2- One can show (see section 3 of this paper) that this is a general 
property of the Liouville boundary structure constants: all they satisfy recurrence identities 
that have the form of difference equations with respect to one of the boundary parameters. In 
particular, using the pentagonal equation satisfied by the boundary three point function , we 
find two difference equations for the latter, one of them quite similar to the equation obtained 
in ^2j. The shifts of the boundary parameters come from the fusion with the lowest degenerate 
boundary Liouville fields. 

The two and three point functions in 2D QG are products of the corresponding Liouville 
and matter structure constants, and therefore have the same dependence on the boundary 



(jlb = fJ-B cosh T. 
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parameters. As a consequence, they should satisfy difference equations of the same form. 

In this paper we will concentrate in the simplest realization of the matter field is as a 
gaussian field with background charge. Then the matter three-point function is trivial (no 
matter screening charges) and the difference equation obtained in Liouville theory should hold 
also for the three-point boundary correlator in 2D quantum gravity. 

In order to check this result experimentally, we need an adequate lattice realization of the 
gaussian field coupled to Liouville. This realization is given by the SOS model on a random 
surface and the related statistical models: the 0(n) JH] and the ADE ^JJEIEI] models of 
quantum gravity. In all these models the matter fields are described as a loop gas on the 
random surface. 

1.2 The discrete approach 

The loop gas realization of 2D QG has been particularly useful in studying problems involving 
two-dimensional geometrical critical phenomena, typically associated with non-unitary CFT 
(see, for example, the review 1191). After being transposed to a random surface, many of those 
problems become easily solvable. The solutions then can be used to obtain information about 
the same problems before coupling to gravity. In particular, the scaling dimensions in the theory 
on the plane are related to the gravitational dimensions by KPZ scaling formula [HIE!- The KPZ 
rule was checked on various microscopic realizations of the 2D quantum gravity in terms of 
statistical models on random lattices and in the last years have been applied successfully to 
evaluate some difficult "flat" geometrical critical exponents associated with percolating clusters 
and polymers [T§] . 

In the loop gas model, the matter fields are represented by non-intersecting loops (which can 
be also interpreted as domain walls) on a randomly triangulated surface. The corresponding CFT 
is non-unitary, with continuous spectrum of the central charge — oo < c < 1 20]. The matter 
correlation functions describe the critical behavior of networks of linear polymers |U 5 . Bulk 
correlation functions of linear polymers have been first calculated in [2Sj using a factorization 
property of the measure over random surfaces. The same method is readily generalized to the 
case of boundary correlation functions. More complicated examples of boundary correlators 
have been calculated in [2U, including the 2-point function for the boundary twist operators 
intertwining between the free and fixed boundary conditions. 

The exact results presented in |12| inspired the work j2S]) where it was shown that the 
boundary two-point function of star polymers on a random disc satisfies, when appropriately 
normalized, the same difference equation as the Liouville boundary reflection amplitude for 
degenerate fields of type (r, 1) or (1, s). It was also noticed that the two-point function of the 
boundary twist operator calculated previously in [23j agrees with the general formula presented 
in ^5] with the corresponding identification of the twist operator as a boundary conformal field. 

Here we apply the approach of |25j to derive difference equations for the boundary three- 
point function and the bulk-boundary correlation function. Comparing the results with those 
obtained in Liouville theory we find that there is perfect agreement between the continuum and 
the discrete approaches. 

The layout of this paper is as follows. In sect. 2 we sketch the world sheet description of 2D 
QG as a theory of a gaussian matter field x with background charge eo, and a Liouville field 4> 
with background charge Q and exponential interaction. We introduce the KPZ primary fields in 
the bulk and on the boundary. Sect. 3 contains a review of the boundary Liouville problem and 

5 This geometrical description of CFT was recently used in the analysis of stochastic conformal maps |22|. 
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the recent results [TSlnailUj . There we summarize the finite-difference equations satisfied by the 
boundary structure constants. In sect. 4 we interpret the difference equations obtained in pure 
Liouville theory from the point of view of 2D QG with gausian matter field. In sect. 5 we give 
the microscopic definition of the 2D QG as a non-restricted height model on the dynamically 
triangulated disc. First we define the height model on a given triangulation and establish the 
equivalence with a gas of loops. We give the loop- gas formulation of the Dirichlet and Neumann 
boundary conditions as well as that of the electric, magnetic and twist operator of the gaussian 
field. In sect. 6 we derive the difference equations for the basic boundary correlation functions 
in 2D QG. First we remind how the loop equation for the disc amplitude can be solved by being 
transformed into finite-difference equation. Then, using the knowledge of the disc amplitude, 
we derive another type of functional identities using a factorization property of the functional 
measure of 2D QG. These identities are most easily obtained for the correlation functions in- 
volving degenerate matter operators with Dirichlet boundary condition. These operators are 
described geometrically as sources of open lines starting at the boundary. We also consider the 
case where both Dirichlet and Neumann boundary conditions are present. Here we derive a dif- 
ference equation for the correlator of two intertwiners between Dirichlet and Neumann boundary 
conditions, which we call boundary twist operators, and a vertex operator with arbitrary charge. 
Our conclusions are presented in sect. 7. 



2 The world-sheet description of 2D QG 
2.1 Liouville Quantum Gravity on a disc. 

Let A4 be the world-sheet manifold, which we take with the topology of a disc. We denote by 
R and K correspondingly the local gaussian curvature in the bulk and the geodesic curvature 
along the boundary 6 . In the conformal gauge g a b( x ) ~ ^ab the fluctuations of the metric are 
described by a local field 4> (the Liouville field) with exponential interaction, which couples to 
the matter field through the conformal anomaly pQ. The conformal gauge necessarily introduces 
the Faddeev-Popov reparametrization ghosts with central charge —26. We will assume that the 
matter is realized as a scalar field \ with background charge eo and central charge c = 1 — 6eQ. 
Such a field is described by an effective action 7 



A[ X , <t>] = j (^t(W) 2 + (VX) 2 + (Q<t> - ie oX )R] + ^ 2b ^j 
M 

+ f C^(Q4> ~ieox)K + fi B e b ^j + ghosts (2.2) 
dM 

where [i and [Ib are correspondingly the bulk and the boundary cosmological constants and the 
background charges are expressed in terms of the Liouville coupling constant b as 

Q = \ + b, eo = \~b. (2.3) 

With the choice Q2.HJI the two background charges satisfy Q 2 — eg = 4, which is equivalent to the 
vanishing of the central charge c to t = + c x + c g h os ts = (1 + 6Q 2 ) + (1 — 6eg) — 26 = 0. The bulk 

6 We normalize the two curvatures so that J M + 2 J dM K = A-k. 

7 Here we choose the boundary interaction that supports the quasiclassical description. The other possible 
interaction is fiB e^ 6 . 
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and boundary cosmological constants are coupled correspondingly to the area A = J d 2 z e 2b ^ z ' s ' 
and the boundary length £ = J dx e b ^ x ^ . 

After conformally mapping the disc M to the upper half plane, the curvature term disappears 

oo 

AM - j *. + (V X ?l + ^) + j dx fis * + g ho St s ,2.4) 

lmz>0 — oo 

and the background charges are introduced through the asymptotics of the fields at spatial 
infinity 

4>(z, z) ~ -Q log \z\ 2 , x(z, z) ~ -e log \z\ 2 . 

The partition function on the disc is defined as the functional integral over the fields x an d <fi 
with action 1)2.2(1 

*(m,Mb) = j\dx\m e- A ^. (2.5) 

To define the functional measure we have to fix the boundary conditions. In order to have 
correspondence with the sum over triangulated surfaces one should impose non-homogeneous 
Neumann boundary condition for the Liouville field 8 

i{d - B)4> = Attub e b ^ x \ 

while the matter field x can satisfy either Dirichlet or Neumann boundary conditions. 

2.2 Bulk and boundary matter fields in the Coulomb gas description 

Bulk matter fields: There is a continuum of holomorphic bulk primary fields realized either 
as electric or as magnetic operators. The electric, or vertex, operators 

V e (z,z) = e 2ie ^ z ~ z) 

are primaries with respect to the energy momentum tensor 

T(z) = -(d X ) 2 - ie d 2 X , f{z) -(B X ) 2 - ie B 2 X 

and have conformal weights 

A e = A e = e(e - e ) = (e - e /2) 2 - eg/4. (2.6) 

The magnetic, or vortex, operators V m (z,z) are associated with discontinuities of the field x- 
The vortex operator with magnetic charge m describes a discontinuity 5x = 4-7rm along a line 
starting at the point (z, z) and has left and right conformal dimensions 

A m = m(m- e ), A m = m(m + e ). (2.7) 

In presence of a background charge the purely magnetic operators have spin. On the other 
hand, only the spinless fields survive in 2D QG, since the local rotations can be compensated 



s Very recently, the Dirichlet boundary condition for the Liouville field studied in I2(i| was also found to have 
interpretation in terms of a sum over random surfaces |27j . It describes the non-perturbative effects related to 
D-instantons and D-particles. The role of the Dirichlet Liouville boundary condition in two-dimensional string 
theory has been previously discussed in |28|. 



6 



by coordinate transformations. Therefore we are led to consider the mixed electric-magnetic 
operator O e ^ m , which have left and right conformal dimensions 

A eim = (e + m)(e + m - e ) A e , m = (e - m)(e - m - e ). (2.8) 

The condition that the spin A — A of such an operator vanishes fixes the electric charge to 
e = eo/2 or the magnetic charge to m = 0. Thus the spinless operators carrying magnetic 
charge are of the form O eQ /2,mi with conformal dimensions 

A eo / 2 ,m = Ae /2,m = TT? - Cq/ '4. (2.9) 

An important set among the primaries are the fields degenerate with respect to the conformal 
symmetry algebra |U| with scaling dimensions 

A rs = {r/b -f- e ° (r, s & N). (2.10) 
These fields can be constructed either as vertex operators with charges 

e rs = ^(e -r/b + sb) (2.11) 
or alternatively as spinless operators with electric charge eo/2 and magnetic charge 

m rs = -(-r/b + sb). (2.12) 

The reflection symmetries e — ► eo — e and m — ► — m act as {r, s} — > {— r, — s}. Nontrivial 
correlation functions of these operators can be constructed by adding a finite number of screening 
operators V\j\, and V-b with conformal dimension one. 

We will be interested only in the magnetic charges of the form m = Lb/2, L £ Z, which have 
a nice geometrical meaning in the microscopic approach. This spectrum of magnetic operators is 
compatible with the assumption that the gaussian field is compactified at radius b: x ~ X + 2irb. 



Dirichlet and Neumann boundary conditions: For the matter field, we will consider 
Neumann (iV) boundary condition 

(8 -B)x = (Neumann) (2.13) 

or Dirichlet (D) boundary condition 

(8 + B)x = (Dirichlet). (2.14) 

These are the only boundary conditions that can be imposed on a gaussian field. 

In order to define completely the D and N boundary conditions, we have to fix the global 
modes. In the case of Dirichlet boundary condition, eqn. 12. 141 means that the field has constant 
value of x( x ) = Xb along the boundary. Thus there is a continuum of Dirichlet boundary 
conditions labeled by the value of xb- Similarly, in the case of Neumann boundary condition we 
should specify the increment of x along the boundary: 5x = / x( x )dx. The Neumann boundary 
condition is related to the Dirichlet one by duality transformation 

b<->l/b, X = x{z)+X{z) ^ X = x{z)-X{z)- (2-15) 
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Boundary matter fields: Boundary fields interpolate between different boundary conditions. 
To completely characterize a boundary operator, one should specify the boundary conditions on 
the left and on the right [22]. The boundary primary fields can be constructed out of vertex and 
vortex operators. There are three types of boundary operators for the gaussian field, N/N, D/D 
and N/D. 

• The N/N boundary fields are the boundary vertex operators 

B e (x) = e iex{x \ 

where x £ R is the coordinate on the boundary. The boundary scaling dimension of this field is 
given by the same formula as for the corresponding bulk operator: 

Af = e(e - e ) = (e - e /2) 2 - eg/4. (2.16) 

• The D/D boundary fields are boundary magnetic operators B m , which introduce a dis- 
continuity 6x = 27rm of the boundary value of the field. The magnetic operator changes the 
boundary value of the matter field from xb to xb + 2irmn. The dimension of such operator is 

A^ = m(m - e ) = (m - e /2) 2 - eg/4. (2.17) 

Unlike the spinless magnetic operator in the bulk, the boundary magnetic operator does not carry 
electric charge. The degenerate electric and magnetic boundary operators have the spectrum of 
charges 

e rs = -(e - r/b + sb), m rs = -(e - r/b + sb). (2.18) 

We will be concerned only by the magnetic charges of the form m = L/2 with L integer, which 
can be realized microscopically as L domain lines ending at the boundary (see sect. 5). 

• Finally, there are the D/N and N/D boundary operators. We call the operators of lowest 
dimension intertwining between Dirichlet and Neumann boundary conditions boundary twist 
operators. For b = 1 (c = 1) the boundary twist operator is well known |30|I31| and its dimension 
is At = 1/16. The boundary twist operator for the Coulomb gas with b < 1 will be introduced 
in sect. 5.4. Its dimension is that of a "half-degenerate" field with {r,s} = {^,0}. 

2.3 Gravitationally dressed bulk and boundary fields 

Bulk operators: We will restrict ourselves to the vertex operators since the magnetic opera- 
tors are vertex operators for the dual field x- I n 2D QG the vertex operators of the matter field 
are dressed by Liouville vertex operators 

V e>a (z, z) = e 2iex(z ~ z) e 2a(t>(z ' 2) (2.19) 

in such a way that the composite operator is a density, i.e. its conformal dimensions are 
A=A=1| 

e(e - e ) + a{Q - a) = 1. (2.20) 
In the language of the 2D string theory, the shifted matter and Liouville charges 

p = e Q - 2e, E = i(Q - 2a) (2.21) 
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are the components of the momentum in the two-dimensional Euclidean target space 
(X®,X l ) = ((f), x)- The balance of dimensions Q2.2UJI then takes the form of the mass-shell 
condition 9 

p 2 + E 2 = 0. (2.22) 

The states with real Liouville energy E appear as normalizable intermediate states, while the 
non-normalizable on-shell fields (|2.2^ describe massless particles in the target space, the closed 
string 'tachyons', which correspond to local operators on the world sheet. The two solutions 
E = =piP of the mass shell condition (|2,22|) describe closed string tachyons with left and right 
chirality 10 

V ^ ) = ^ r(l^S) Ve ' a± > e = \(e -P), a ± = l -{Q T P). (2.23) 

In this way, there are two possible gravitational dressings of the matter field vertex operator 
with charge e. For instance, for the identity operator e = we have two possible Liouville 
dressings: the Liouville interaction a = b and dual one a = 1/6. The solution Vp with the 
smallest Liouville exponent describes physical a physical state, 

(P>0) (2.24) 
(P < 0) 

while the larger exponent is usually referred as 'wrongly dressed' state, 

IV - r {', , (P > 0) (2.25) 

^ (P<0). v ; 

The condition a < Q/2 for the physical fields is also known as Seiberg bound |32j . The wrongly 
dressed states are not realized in Euclidean 2D QG, but they will be important when comparing 
the microscopic quantum gravity with Liouville theory. The correlation functions of the oper- 
ators Vp are obtained from the correlation functions of the 'physical' operators Vp by analytic 
continuation with respect to P beyond the point P = 0. 

The gravitational dressing of the identity operator (e = or P = eo) is the Liouville inter- 
action V eo = e 2 ^; this implies that we have chosen b < 1 in the action l2~2l 





Boundary operators: Similarly to the bulk fields, the boundary fields in 2D QG are con- 
structed as matter and Liouville vertex operators 

B et f)(x) =e ,exW+Wa; ». (2.26) 

In the string-theoretical context, the vertex operators with complex momenta 

e = I eo -P, P = ±Q-iE (2.27) 

9 Here we assume that the Liouville field describes the time direction of the target space of the Euclidean 2D 
string theory. 

10 Since we will be eventually interested in the comparison with the discretized theory, we choose a normalization 
which eliminates the "leg factors" in the correlation functions |3I4| . Below we will normalize the boundary fields 
following the same principle. 
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appear as normalizable intermediate states of the open string amplitudes while the on-shell fields 



are the open string tachyons. The physical boundary operators 

(P>0) 
(P<0) 

are related to the fields with the 'wrong' dressing 





(P>0) 
(P < 0) 



by analytic continuation through the point P = 0. 

We will denote the degenerate matter fields dressed by gravity as 

B rs = B Pra , Prs=r^ h -s h -. (2.31) 

As we will see later, the degenerate fields £>l+i,i and B\.l+i play special role in the microscopic 
approach. 

In order to define completely a boundary operator, we have to specify the left and right 
boundary conditions both for the matter and Liouville fields. The Liouville boundary conditions 
are labeled by the values of the cosmological constant on both sides. 

Let us remark that here both dressings can have physical meaning. In particular, the bound- 
ary Liouville interaction can be e b ^ as well as e^' b , the two exponents corresponding to the 
two possible dressings of the identity operator. In the Liouville literature only the first case 
(which was chosen as the boundary term in the action I2.2j) has been considered. In this case 
the dimension of the boundary is one, i.e. half the dimension of the bulk: ~ /W 2 and the 
boundary can be treated quasiclassically. In the second case /j,b ~ ^ x l 2h2 , which corresponds 
to the "wrong" dressing of the identity boundary operator, the boundary has anomalous fractal 
dimension larger than one. It has been shown that both boundary Liouville interactions can be 
realized microscopically |16| ll8j. This is why it is useful to introduce a new scaling exponent v 
having the meaning of inverse fractal dimension of the boundary 



~ H (2-32) 



where 



v = 1 for boundary interaction e ^ 
v = b 2 for boundary interaction e^ 6 . 

The two possible values of v are dual to each other in the sense that if the classical boundary 
{y = 1) corresponds to a Dirichlet boundary condition, then the fractal boundary (y = b 2 ) 
corresponds to Neumann boundary condition for the matter field x> an d wi> ce versa |24j . 
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Correlation functions and KPZ scaling: Consider the correlation function of n boundary 
and m bulk fields as 11 

W^:;^ m \p B1 ,...^ Bn ^) = (Bpf 1 ,Ms2 . . .Bp^ n '^ Bl V^.-.V^)^, (2.33) 

where ( ...) denotes the functional integral with respect to the matter and Liouville fields. By 
the conformal invariance this correlation function has the form 

^F^^- 1 ^,...,^- 1 ^), (2.34) 
where the exponent 7 depends linearly on the target space momenta 

1 m 1 n 

± \ "V , bulk. ± \ "V , bound. , , 

7 = l-- 7str -^(l-^. )--J2(l-S Pk ). (2.35) 

3=1 k=i 

Here u\, u n are the exponents (j2.33j) for the n segments of the boundary, 

7str = 1 - 1/b 2 < (2.36) 

is the string susceptibility exponent, 

bulk \K\ 1 bound \P\ 1 

ok = + 2 7str 6p = ~b~ + 2 7str ^ ' 

are correspondingly the gravitational dimensions of the bulk and boundary fields, and F is a 
scaling function. The gravitational dimensions 5 are related to the flat dimensions A and the 
conformal anomaly c of the matter field by the KPZ formula 

A = 6(6 ~ 7str) , c = l-6 7L ■ (2.38) 

1 - 7str 1 - Tstr 

Using the results obtained in boundary Liouville theory and reviewed in the next section, we 
can evaluate instantly the boundary correlation functions in 2D QG that involve three or less 
fields. In this case there is no moduli integration and the correlation functions of the matter 
and Liouville fields factorize and the correlation function in 2D QG is equal (up to a numerical 
factor) to the corresponding Liouville structure constant. There is a subtlety related to the fact 
that only physical fields are present in 2D QG. It is resolved by using the reflection symmetry 
to express the correlation functions of wrongly dressed fields (whenever they appear) in terms 
of the correlation functions of physical fields. 



3 Correlation functions in boundary Liouville theory 
3.1 Review of boundary Liouville theory 

In this subsection we present a brief review of some of the exact results for Liouville on a disc, 
obtained recently in |lllll2j . In the upper half plane the Liouville theory is defined by the action 




UHP R 

with asymptotics at infinity (j)(z, z) ~ —Q \og{zz) and Neumann boundary condition on the real 
axis i(d — d)4> = 4717x5 e 6( ^ '. The theory is conformal invariant with central charge cl = 1 + QQ 2 . 

To completely define the correlation function, we should also specify the matter boundary conditions, which 
is understood below. 
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Bulk fields: The fields V a (z, z) = e 2a ^ z ' z ^ are primaries with respect to the energy momentum 
tensor 

T(z) = -(d<j)) 2 + Qd 2 cp , 

f(z) = -{3(f)) 2 + Qd 2 4> 

and have conformal weight A a = A a = a(Q — a). Because of the invariance a — > Q — a, one 
identifies the operator V a with its reflected image Vq_ q : 

V a (z,z) = S(a)V Q - a (z,z) (3.2) 

where we used the bulk reflection amplitude 

(7r^(fe 2 ))(q- 2 ")/ 6 7 (2a6 - b 2 ) 
S{a) = b 2 7 (2 -2a/b + l/b 2 Y (3 ' 3) 

and 7(1) = r(x)/T(l — x). The reflection amplitude satisfies the unitarity condition: 

S(a)S(Q -a) = 1. (3.4) 

An important set among the primaries are the fields V_ n \ ) j2, n G N, which are degenerate with 
respect to the conformal symmetry algebra and satisfy linear differential equations 9 . For 
example, the first non trivial case consists of a = —6/2, and the corresponding operator satisfies 



8 2 + T(z) V_ b/2 = 0, (3.5) 



as well as the complex conjugate equation. It follows from these equations that when one 
performs the operator product expansion of one of these degenerate operators with a generic 
operator, then the OPE is truncated 9 . For example: 

V- b /2V a = c+V a _ b/2 + C-V a+b/2 ■ (3.6) 

The structure constants c± are special cases of the bulk three point function, and can be com- 
puted pertubatively as Coulomb gas (or screening) integrals [SS1E1]- One can take c + = 1, as 
in this case there is no need of insertion of interaction, whereas c_ requires one insertion of the 
Liouville interaction — [i J e 2b( ^d 2 z, and 

c- = -li /^(^(0)F_ 6/2 (l)e 26 ^)F Q _ a _ 6/2 (oo) 



= -nj d 2 z\z\- 4ba \l-z\ 2b2 

7T7(26a - 1 - b 2 ) 
_/i 7(-6 2 )7(26a) ' 

In the first line, we used the property of invariance under global transformations to set z\ = 
0, Z3 = 1, £4 = 00, and in the second line < cft(x)(j)(y) >= — log \x — y\. 

Similarly, there exists also a dual series of degenerate operators V- m /2b with the same properties. 
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Boundary fields: There is a one-parameter family of conformally invariant boundary con- 
ditions characterized by different values of the boundary cosmological constant For later 
convenience we will parametrize \xb m terms of a dimensionless angular variable a defined as 
follows Q2]: 

cos (2nb(a - Q) ) = ^ ^(vrft 2 ). (3.7) 

The meaning of this parametrization is that it unfolds the branch point of the observables at 
Hb = —^°b = v/V sin7r6 2 . We shall denote the boundary operators as BZ 2(T1 (x); they have 
conformal weight = (3[Q — f3) and are labeled by two left and right boundary conditions ai 
and o"2 related to fis 1 and /ie 2 by (|3.7() . Again, because of the invariance /3 — ► Q — /3, one should 
identify the field B^ 2CT1 {x) with its reflected image B^^(x): 

B^(x) = D(J3,<t 2 ,<ti) B^(x) (3.8) 

The boundary reflection amplitude D{(3, a 2 , ci) satisfies the unitary condition 

D((3, a 2 , ai)D(Q - (3, a 2 ,a 1 ) = 1 (3.9) 

and gives the nontrivial piece of the boundary two-point correlation function, see the next 
section. 

As it was argued in JI2|. the degenerate boundary fields B^ b ^ 2 , n £ N, have truncated 



operator product expansion with all primary fields if the left and right boundary parameters are 
related by a 2 — o\ = — — ( n ~~ 2 ) 6 ; , , For example, for a 2 = a\ ± | 

Bf^B^ = c%B;t% 2 + c ± B^ b ) 2 , a 2 = a 1 ±b/2 (3.10) 

where, as in the bulk situation, c± are structure constants and are obtained as certain screening 
integrals. In the first term of (|3.10|) there is no need of screening insertion and therefore ci can 
be set to 1, whereas the computation of requires one insertion of the boundary interaction 
—\ib f e b ^ x ^dx and was explicitly computed in ^2] with the following result 

c - = 2 ( /^ToT r(26/3 - b 2 - l)r(l - 26/3) x 

sinirb(P ± (<7i - <r 3 ) - 6/2) sin tt&(/3 ± (<ti + a 3 - Q) - 6/2). (3.11) 



Similar properties have the fields B_^ n J b with a 2 — a\ = — ^ 



n n— 2 ra 
26 ' ' ' 26" 



Self-duality: The observables are invariant with respect to the duality transformation 6 — ► 1/6 
provided the dual cosmological constant p, is related to fx as 

vr / i 7 (l/6 2 ) = (vr M7 (6 2 )) 1 ^, (3.12) 
and the dual boundary cosmological constant is defined as follows 



^(t^-I^ vlv sin 5' (3 ' 13) 
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3.2 Correlation functions 



In order to characterize LFT on the upper half plane, one needs to know the following structure 
constants: 

1. bulk one-point function 

\z — z\ 

2. boundary two-point function 

S(P 2 + Pi-Q) + D(/3 u a 2 , aJSifa - Pi) 



3. bulk-boundary two-point function 12 

R a (a,f3) 



\Z — Z\ p \Z — X\ ^ 



(V a (z,z)B^(x)) 
4. boundary three-point function 

B ' 7 ^ :> (xi)B a ^ 1 (x2)B (T ^ 2 (x 3 )j = Ai+A2 _a 3 |A2+A 3 -Ai I 1A3+A1-A2 

' F2l| \X32\ F3l| 

By its definition, this function is invariant under cyclic permutations. In the following we 
will also use the notation 

^y(<T2CT3(Ti)/33 ^(0"20"30"l) 

°/3i/3 2 - U Pi/3 2 ,Q-l3 3 - 

All correlation functions considered above are invariant w.r.t. the duality transformation 
b^ 1/6. 



3.3 Explicit expressions and functional relations 

Below we list the explicit representation for the structure constants, and the functional equa- 
tions they satisfy. Except for the boundary three point function, they have been computed by 
considering an auxiliary correlation function containing a degenerate bulk or boundary field: 
this leads to solvable functional relations for the correlation function sought for. This trick was 
first introduced in [33] , and used extensively in jHHl^JETlE] • For each of the equations listed 
below, there exists a dual equation, obtained by replacing b by 1/6. It is important to note that 
the functional relations are of two types: one type involves shifts on boundary parameters and 
momentum, the other type involves shifts of boundary parameters only. 

1. bulk one-point function 
Expression MMSty : 

2 9 (Q-2a) 

Ua{0i) = -(TTflj(b 2 )) 2b X 

xr(26a - 6 2 )r(^ -77-1) cos[^(2a - Q)(2a - Q)} (3.14) 
b b z 

12 The bulk one point function is a special case of the bulk-boundary coefficient with f) = 0. 
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Reflection property: 

U a {a) = S{a)U a {Q - a), 
where S(a) is the bulk reflection amplitude 1)3. 3 j) . 

Functional equation: 

U _ k (a) + U , k (a) = 2 cos irb(2a - Q)UJa). (3.15) 

2. boundary two-point function 
Expression )1^ : 



^ F b (2(3 - Q) S b (a 2 + Q-1- (3)S b (2Q - - at - a 2 ) 
X T b (Q-2(3) s b (P + *2-<ri)S b (0 + a 1 -* 2 ) 

Reflection property: it takes here the form of the unitarity condition 

£>(A<7 2 ,ai) D(Q-P,o- 2 ,o-i) = l. (3.17) 

Functional equations: The first of these equations can be seen as a consequence from the 
different possible factorizations for an auxiliary boundary three point function containing 
a degenerate boundary field < -Bg+6/ 2 (x3)-B-6/ 2 (;c 2 )-Bg(xi) >, see equation (4.3) of [T2] . 



D(/3, o-^o-! - 6/2) - D((3, a 2 , a x + 6/2) 
-2tiT(1 - 26/?i) 



I 1 



T(2 + 6 2 - 26/3i) 

One can also find 



vt7(— 6 2 ) 

sin7r6(2cri - Q)D(/3 + 6/2, a 2 , a x ). (3.18) 



sin 2 7r6(<7i - o 2 - [3 + ^)D([3, a 2 + ^, cii) 

- sin 2 7r6(ai + ct 2 - /3 - ^)£>(/3, cr 2 - -, a x ) 

= sin7r6(2o-i - Q) sin 7t6(2ct 2 - Q)D(/3,a 2 ,o-i + ^), (3.19) 

and when the two boundary parameters are equal o\ = <7 2 = a: 

sinvr6(2a - Q - (3)D a _ b/2 (p) + sinvr6(2a -Q + 0)D a+b/2 ((3) 
= 2sinvr6(2CT - Q) cos vr6(/3 - Q)D a (/3). 

3. bulk-boundary two-point function 
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Expression f 14[\3 



R a (a,P) = 2ir(ir^(b 2 )b 2 - 2b2 )T b (Q-^-P) 

TUQ - (3)T b {2a - f3)T b (2Q -2a- [5) 
r b (Q)T b {Q - 2p)T b (/3)T b (2a)F b {Q - 2a) 

loo 

dp e 



2^(2.-0) S b ( P + (3/2 + a- Q/2)S b (p + (3/2-a + Q/2) 



S b {p -(3/2-a + 3Q/2)S b (p -(3/2 + a + Q/2) ' 

Reflection properties: 

R a (a,f3) = S(a)R a (Q-a,f3), R a (a, (3) = D(p,a,a)R a {a,Q - /3) 13 . 

Functional equations: The two functionals relations written below are nothing but the 
Fourier transform of the equations (21) found in |14j . We recall that these equations arise 
from the different possible factorizations for the auxiliary three point function containing 
a degenerate boundary field < V a (0)Bp(l)B_ b / 2 (r]) >• We find: 



R„-k{a,P) - R„,b(a,/3) = sin7r6(2<T - Q)R a (a,P + b) 

CT 2 ' 2 



X 



n \ 1/2 r(i - 26/3)r(i - b 2 - 2bp) 

-2tt 



TT-f(-b 2 ) J T 2 (l - b(3)T(l - 6/3 - 26a)r(l - b(3 + 2ba - bQ) ' 

(3.20) 

as well as 

sin7r6(2cr - Q - (3)R a _ b/2 (a, (3) +sm-nb(2a -Q + (3)R a+ b(a, (3) 
= 2sin7r6(2cj - Q)cosvr6(2a - Q)R a (a,/3). 

(3.21) 

One notices that the special case f3 = in the latter equation reproduces equation ()3.15|) . 

Particular case: When 2a = (3, one can show the following relation that links the bulk- 
boundary structure constant to the boundary reflection amplitude: 

R a ((3) = ies 2a= pR a (a, (3) 

r 2 (Q-/?) 

171 

4. boundary three-point function 

13 The first of these relations is trivially verified. We checked the second one is also true, as it was not done 
in [H]. 



2tt b ^ ^ D{(3, a, a)re Sx=0 T b (x), (3.22) 
T(((3)T b (Q) 
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Expression llty : 



r b (Q + P2-P1- Ps)r b (Q + fh - fa - 02) 

X r 6 (Q - 20 x )T b {Q) 

T b (2Q -fa-fo- 3 )T b (0 2 + 03 - Pi) 
X T b (20 3 - Q)T b (Q - 202) 

S b (03 +CF1- &3)S b (Q + 03 ~ 03 ~ gl) 
5a(^2 + <7"2 - CT3)S b (Q +02 -OS- (T 2 ) 

ioo 

X i i S S.iVr + s)S b {V 2 + s)S b (V 3 + s)S b (Q + s) 1 ' ' 

—ioo 



and the coefficients Ui, 


Vi 


and i 


= 1, 


' ' ' 1 


4 read 






Ui = 




+ 02 ~ 


-A 




Vi = 


Q + 0-2- 


0-3-0! +0 3 


u 2 = 


Q 


— (71 + <7 2 - 


-01 


v 2 = 


2Q + a 2 


-(73-01-03 


u 3 = 


02 


+ (72 " 


- <73 




v 3 = 


2(7 2 




u 4 = 


Q 


-02 + 0-2- 


~ 03 









Reflection properties: It was shown in that 

^-,(0-20-30-1) _ /-((o"20"3(7i)Q— 3 rtf a „ „ \/^»(^2Cr3(7i)/33 

C 010203 =G /3i/3 2 =V{0 3 ,O-2,O- 1 )Cp i p 2 

as well as two other similar equations obtained by cyclic permutations. 



Functional equations: It is shown in the appendix B that the following two equations are 
a consequence of the pentagonal equation satisfied by the boundary three point function 1 

^,{^2,0-3,^1—2) _ ^(°"2, 0-3,0-1 + !) _ ^-,(02,0-3,01) 

01 ,02,03 U 01,02,03 - ^hfo+kfo+k 

( V \ V2 

x — 2ir I ■; — -rp- sin7r6(2(7i — Q) x 



and 



7T7( — 6 2 ) 

T(l - 2b0 2 ) T(l - 26/3 3 ) 



T(l - b(0 2 + 03 - 0i)) r(2 - + & + 03 - b)) 



(3.24) 



sinvr6((72 - cri - /% + -) sin7r6(<7 3 - (7i - 2 + -) ^SiXX*"^ 

- sin vr6((7 2 + (7i - 3 - -) sin 7r6(<7 3 + cri - /? 2 - ~) C^^/C 2 } 

= sinvr6((72 + a 3 - 0i) mnrb{2a x - Q) ^+*£ s+ *' ai) . (3.25) 



14 We recall that this fact follows from the consistency condition that expresses the associativity of the product 
of four boundary operators |38II39| . 
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These equations reproduce equations (|3.18j) and (|3.19[) in the case 0\ = 0, 02 = 03, 
and 02 = ^3, where the structure constant C^fa''* 1 reduces to the boundary reflection 
amplitude D(0, 02, ci). 



3.4 Important cases 

To compare the Liouville approach with the matrix models approach, it is useful to compare the 
dependence with respect to the boundary parameter a of some important correlation functions, 
such as the one point function of the boundary cosmological operator B? a , the bulk-boundary 
structure constant with one boundary cosmological operator inserted R a (a, b), and the boundary 
three point function with three boundary cosmological operators inserted. 

• < Bl u >: 

It amounts to evaluate the quantity R a (0,b) directly from its integral representation. In 
sending a — > 0, some of the poles of the Sb functions will cross the imaginary axis, so 
one has to deform the contour accordingly. Picking up the residues at poles leads to the 
following expression: 

< Bl a >~ / u^costt&~ 1 (2o--Q). (3.26) 

In other words, the one point function of the boundary cosmological operator B^ a is 
proportional to the dual boundary cosmological constant \xb- 

• R a (a, b),a 7^ : 

Setting = in (EHffll) provid es a relation between R a (a,b) and the one point function 
Ra(a, 0) = U a {a). Tor the time being, we are only interested in the dependence of i? CT (a, b) 
and U a {a) on <r, which is found to be: 

RM,b) ~ ^(^) ^^(2a- g)(2a-g)] 

sin7ro(2<7 — Q) 



C bbb 



(<720"3°"l) . 

bbb ___ 

It is not straightforward to derive it directly from ()3.23|) . but one can check that the 
following expression 



is a solution of the functional equation (|3.25|) . 



4 Boundary correlation functions in Liouville quantum gravity 

In order to carry out the comparison with the discrete approach, we need to take into account 
the matter field and also normalize the fields and the coupling constants in a way to match with 
their microscopic realizations. Among the results listed in the previous section it is sufficient 
to consider eqn. (|3.24|) together with its dual since it contains all the information we need 
about the boundary correlators. We will show that eqn. (|3.24[) together with the duality and 
reflection properties is sufficient to determine, up to a normalization constant, the one-, two- 
and three-point functions in 2D QG. 
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Since we are considering the simplest case of a free field with no screening charges, the role 
of the matter field is reduced to imposing the neutrality condition for the charges of the vertex 
operators. We assume Neumann boundary conditions for the matter field along all segments of 
the boundary. Thus the difference equations in Liouville theory hold also for the physical fields 
in 2D QG. 

The only subtle point here is that the equations do not close on the set of the physical fields. 
Repeated application of the difference operator produces eventually a "wrongly dressed" field. 
In this case we will need the boundary reflection amplitude to express the "wrongly dressed" 
field through a physical one. 

Once closed equations for the three-point functions are found, we will formulate them entirely 
in terms of the observables that can be measured in the microscopic theory. For this purpose 
we evaluate the boundary one- and two-point functions in 2D QG. 

4.1 The rescaled bulk and boundary cosmological constants 

We will redefine the bulk and boundary cosmological constants fi and [ib according to the nor- 
malizations (|2.23j) and (|2.29j) for the bulk and boundary operators. The new bulk and boundary 
cosmological constants 15 A and z are defined as the couplings of the operators V$ and 
respectively: 

,, [ n 2b<t> _ a v (+) b<f> _ »(+) 



/x ^(1-62)' m r(i-&2)" (429) 



This gives 



We will also introduce the boundary parameter r of [Toj related to a by 

'-t-'sb- (4 - 30) 

Now the boundary cosmological constant is parametrized as 

z = M cosh r, M = VX. 

4.2 On the duality property of 2D QG 

It follows from (|3.12[) - (|3.13j) that the duality 6 — ► b = 1/6 transforms the new bulk and boundary 
cosmological constants into 

A = A 1 / 6 ' , z = M l l b ' 2 cosh r/b 2 (4.31) 
where the dual boundary cosmological constants are defined by 

A V_ eo = p. e^l\ ~z = ii B e^. 

If we denote 

z = Mcoshf, 

then the duality transformation b = 1/6 is formulated as 

A = A 262 or M = M 1 / b \ f = r/b 2 . (4.32) 
15 Here we use the conventions of ref. |18|. 
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The duality acts on left and right chiral tachyons as 



V 



(±) 



V 



(=F) 
-Pi 



[63] nT2 , 



(4.33) 



where the bar means complex conjugation. In the (non-unitary) realization of 2D QG based on 
the gaussian field the duality 6 — > & = 1/6 is not a symmetry of the Hilbert space of the theory. 
It acts in a pair of CFT characterized by matter background charge: eo and e"o = —&o- 
The duality property of the correlation functions thus reads 



Z r (Ki,...,K m ) ,~ 



W, 



(n,...,f n ) 



W D . o (ri,...,r n ). 



--Pi,. 



4.3 The boundary reflection amplitude in 2D QG 



The open string tachyons of same momenta but opposite chiralities are related by the reflection 
amplitude 



[B 



(+JlTlT2 
P i 



(— )j7lT2 



D p + (t 1 ,t 2 ) [4 +) ] 



(4.34) 



The latter is obtained from the Liouville reflection amplitude (|3.16|) by taking into account the 
leg factors in the definition (|2.29|) . Indeed, the conventions used in sect. 2 are such that the 



reflection amplitude for the matter gaussian field e e 
of the rescaled reflection amplitude is 



D (e -e)4> ■ 



is equal to one. The expression 



1 M 2P/b S b [%+P 



2-Kb 



Sh 



+ p 



bS b (2P + b) a(Q P .- 71+7* 

JJ, l o 1 2-Kb 



^ I 2 

oo 



s b % 



1 fdy 



b M T b (2P + b) eXP \2j y 

o 



sinh(7rPy) cos ^P- cos 



yr-2 



sinh(7ry/6) sinh(7r6y) 




— • (4-35) 



The amplitude Dl 



(D p + ) 1 has the symmetries 



DpDZp 



1 sin^P 
¥ sin 2irbP ' 



D p + DZp- 



, sin 2ixbP 
sin^P 



and 



(4.36) 



where the function Dp~(j\, r 2 ) is defined by (|4.35|) with b replaced by b = 1/6, and fi |2 
Remarkably, the function (|4.35|) satisfies difference equation (|3.18|) and its dual without the T- 
functions on the r.h.s.: 



(4.37) 

n, 2 /b 2 . 



a 



sin ( irb 2 — ) Dp~ (ri,r 2 ) = Msinhn (ti,t 2 



(4.38) 



sin 7T 



d_ 

dr\ 



D 



--(ti,t 2 ) = M 1 '* S inh (p) D-+^( n ,n) 



(4.39) 



Note that the reflection amplitudes for the degenerate momenta P mn = \{m/b — nb) are 
given by rational functions of z = Mcoshr and z = M 1//fe2 cosh(r/6 2 ). They can be obtained 
from the reflection amplitudes for P = 



bS(b) 



Di 



1 



S(l/b) 



1 

Y 2 



20 



by applying m times eq. 1)4.39)1 and n times eq. (|4,38[) . Let us calculate the reflection amplitude 
for the Liouville interaction, which will play important role in the sequel. We have 

D b/2 = ■ i ui\ ( Zl + Z2 ) 

+ _ 6 2 sin(7r6 2 ) z\ - z 2 

v -i/2b+b/2 ~ ~ sin ( 7r /52) - _ Z2 

)+- 

1/26-6/2 

where in the last line we used the property (|4.36|) 



or,,-,,, - ( 44 °) 



4.4 Difference equations for the three-point function 

Consider the boundary three-point function, which factorizes to a product of a matter and 
Liouville components 

Wp^Psin^n) = ( [B Pl Y^[B P2 Y^[Bp 3 Y^ } disc = c™?£c%%%. (4.41) 

In absence of screening operators 16 , the gaussian matter field correlator is non-zero only if the 
three electric charges satisfy e\ + e 2 + e 3 = eo or, in terms of momenta, 

Pi + P2 + P3 = ^e . (4.42) 

In this case the matter structure constant does not depend on the matter boundary conditions 
and is equal to one. Indeed, the difference of the matter boundary parameters adjacent to 
a vertex operator is fixed by the charge of the operator. The overall shift of the boundary 
parameter is just a global mode and does not change the three-point function. 

Let us assume that Pi > 6/2, P 2 < —b/2, in which case Bp 1 = and Bp 2 = B p \ 

Taking into account the normalization ([2.28ft and using the neutrality condition 1)4. 42|) , we write 
eqn. (f^M)) as 

Wp u p 2> p 3 {T 1 ,T2,T 3 +ib 2 7r)-Wp u p 2: p 3 (T 1 ,T 2 ,T 3 -ib 2 TT) = 2iM sinhr 3 Wp 1 _ b / 2t p 2+b /2,p 3 (n,T2,T3). 

(4.43) 

Remarkably, all T-functions disappeared. This allows to write the difference equation in the 
operator form 

Ap 2!Pi (t 3 ) W Plr P 2! P 3 (Tl,T2,T 3 ) = (4.44) 

where Ap 2] p 1 (rs) is a difference operator associated with the segment (12) with boundary 
parameter T3 and momenta P 2 and P 2 at the ends: 

Ap 2>ft (r 3 ) = sin (b 2 8r 3 ) - Msinh(r 3 ) et^-^) ^ > ^, P 2 < (4.45) 

In a similar way, the equation dual to 1)3.24)1 can be reformulated as an equation for the three- 
point function ()4.41j) with Bp x = Bp and Bp 2 = Bp^ under the condition that Pi < —1/26 and 

16 The argument that follow can be of course carried out within the Coulomb gas formalism in its full generality, 
which involves an arbitrary number of screening charges. Then the difference equations will depend also on the 
boundary three point function of the matter boundary fields. 
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P 2 > 1/26. The resulting difference equation has again the form (|4.44|) with difference operator 
given by 

Ap 2)Pl (r 3 ) = sin (d T3 ) - M 1 ^ sin (r 3 /6 2 ) e^"^ (p x < -1 P 2 > I") (4.46) 

The difference operator in (|4.44j) involves translations of the momentum by 6/2 or 1/26, and 
after its repeated application on gets on the r.h.s. a correlation function in which one or two 
operators get wrong Liouville dressing. In order to obtain closed set of equations we should 
express the "wrongly dressed" fields (|2.25|) in terms of the physical ones using the boundary 
reflection amplitude. This gives again an equation of the form (|4.44|) . with a difference operator 
which depends on all three boundary parameters: 

-if Pi € oo] and P 2 G [-^,0]: 

AjW 7 *) = sin(vra T3 ) - M 1 ^ 2 sinh (r 3 /6 2 ) D~+ (t u t 2 ) eM 9 ^^) (4.47) 



ifP lG [0,i] andP 2 G[-i,0]: 

Ap 2iPl (r 3 ) = sin(7ra r3 ) - M 1 ^ sinh (r 3 /6 2 ) P+~ ^n, r 3 )P p + (r 3 , r 2 ) e^-^i). (4.48) 



Pl ~2b P2_r 25 

Similarly we get the equations for the intervals Pi G [— |,0],P2 G [f ,00] and Pi G [0, f],p2 G 
[-1,0]. 

It is possible to express these equations entirely in terms of the observable in the microscopic 
theory. Indeed, the hyperbolic sine in the numerator of the second term is proportional to the 
one-point function (j3,26|) and the reflection amplitudes in the denominator are proportional 
to the boundary two-point function. In order to fix the coefficients, we will need the exact 
expressions of the boundary one- and two-point functions in 2D QG. 

4.5 The boundary two-point function 

The boundary two-point function in 2D QG 

W Pu p 2 (rur 2 ) = ( [B p J^[Bp 2 Y^ ) disc 
can be obtained from the boundary three-point function using the relation 

dWp ^i,r 2 ) = { [Bp ^ [Bi _^ [Bp ^ )d . sc = W PuP ^ eo (r U T 2 ,r 2 ). (4.49) 

As a consequence, it satisfies the difference equation (j4.43|) . (|4.47|) as well as the dual equations 
in the corresponding momentum intervals. By the neutrality condition (|4.42j) the two-point 
function vanishes unless Pi + P 2 is equal to zero 17 . We will denote 

Wp(ri,T2) = W- P (n,T 2 ) = Wp- P (n,T 2 ). (4.50) 



17 In a more general setting when screeneng charges are allowed, this condition will be relaxed to Pi + P2 
i/2b-nb/2. 
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This function is proportional to the boundary reflection amplitude, but the evaluation of the 
exact coefficient is not easy in CFT because of the zero mode associated with the residual 
dilatation symmetry. This is why we will evaluate it indirectly, by solving the difference equation. 

The function Wp{j\, t 2 ) should satisfiy (|4.44l) for any P, with the corresponding identification 
of the difference operator. This gives two independent equations for each value of P > 0. For 
sufficiently large P the function Wp(t\,t 2 ) satisfies 



sin 



(b 2 d/dn) - Msin(n) e ~f a / 9p ) W p {t x ,t 2 ) = (\P\ > -6/2) (4.51) 



and 

r sin(d/dn) - M 1 ^ 2 sin(n/6 2 ) e-^ d ' dP ) W P (n,T 2 ) = (\P\ > 1/26). (4.52) 



These are exactly the equations (|4.39|) and ()4,3<S|) for the boundary reflection amplitudes D H 
and D '". Therefore Wp should be proportional to both amplitudes, with P-dependent coeffi- 
cients. The common solution is determined up to a multiplicative constant k = k(6): 

Wp(Tl ' T2) = Kb \ s ^P/b)\ D \ p \ {TuT2) = Kb \sHilpb)\ Dz h {T ^- < 4 - 53 ) 

The duality transformation acts, according to ()4.37j) . as 

Wp(n,f 2 ) = —Wp(T l ,T 2 ). (4.54) 

For \P\ < 1/26 or \P\ < 6/2 the two-point function satisfies non-linear difference equations, 
which follow from (jQ7jl and (|4"i%|) : 



2 



M 1 ^ 2 sinh (n/6 2 ) / 1 



sinfiTroLWi , P] (n , r 2 ) = — .„ p . . — V "V " ' \\P\<ir\ (4-55) 



and 

sin(z7r6 2 a T1 )m_ |p| (Ti,r 2 )= , . - — ^ . „ p . , "VT" v V ( \P\ < ~ ) ■ ( J.oG) 



2 



«g Msinh(ri) / 6 

1 p| 1 ' ' ' ' ' |sin(27r6P) sin ( *f) \ W P (n , r 2 ) V 2 



4.6 The boundary one-point function (the disc loop amplitude) 

The boundary one-point function, or the disc loop amplitude, is defined as 

W(t) = ( [B eo ] TT ) disc = - — , z = Mcoshr 

where $ is the partition function of the disc with fixed boundary condition. It is equal, up to a 
constant factor, to the Liouville amplitude (|3.2fij) . 

As in the case if the two-point function, the direct calculation of W(r) from CFT is difficult 
because of the zero modes associated with the residual global conformal symmetry. We can 
alternatively evaluate W(t) using its relation with the boundary two-point function 18 



18 The easiest way to prove that is to perform a Laplace transformation and express both sides in terms of the 
disc amplitude with fixed length £ = J dx^^ . 
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The explicit expression for W^^Ti, r 2 ) follows from 1)4.40(1 : 

M 1 /" 2 - 1 cosMn/^-cosMVfr 2 ) 

W eo/2 ri,r 2 = . . 4.58 

' |sm(7r/o z )| cosh ri — cosh T2 

Comparing ((4.58(1 and 1(4.57(1 we get the expression for the disc loop amplitude 

\ l Ml/b2 cosh(r/6 2 ) 

W(r) = 6k & — /y ; . 4.59 

| sm7r/o z | 

In this section we are considering pure Neumann boundary conditions for the matter held. 
The loop amplitude is proportional to the volume of the momentum space, which is 2b. If we 
normalize the amplitude by the volume of the momentum space, the factor b in the denominator 
will disappear. Anticipating the result of the comparison with the discrete theory we will fix 
the value of the constant k: 

k=\. (4.60) 
b 

Note that the normalized amplitude is not self-dual. The T-dual amplitude is given by 

WW (4.61) 
tr |sm7r6 z | tr |sin.7r£r| 

4.7 The difference equations in terms of loop observables 

To compare with the microscopic approach we should translate the finite-difference equations in 
terms of boundary observables. 



a) when Pi G [35,00] and P 2 G [—w, — 26J • 



00, 



Wp 1i p 2j p 3 (tx,t 2 ,t 3 + in) - W Pu p 2 ,p 3 (t 1 ,t 2 ,t 3 - in) = 

[W(t 3 + in) - W(t 3 - in)} W Pi __± P2+± _ p (n,T 2 ,T 3 ) (4.62) 

al) when P x G [35,00] and P 2 G [-35,0]: 

Wp l! p 2 ,p 3 (ri,T2,T 3 + in) - W Pu p 2t p 3 (T 1 ,T 2 , t 3 - in) = 
1 W(t 3 + in) - W{t 3 - in) 

1 . 7^ n7r ,„ A, 777 7 ; W P 1 P+ j_ P jTi,r2,r 3 ) (4.63) 

|sin(^ - 2f |P 2 |)| W^_|p 2 |(T3,Ti) 1 2i >' 2+ 2i " 3 

a2) when P x G [0, 35] and P 2 G [-00,-35]: 

^A.ft.ftCn, T2, r 3 + wr) - Wp u p 2! p 3 (t 1 ,t 2 ,t 3 - in) = 
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W(t 3 + 27r) — W(t 3 — in) 



|sin(2vr6Pi 



1^1(7-3,72) 



w, 



(4.64) 



a3) when P 1 e [0, 4] and P 2 G [-4,0]: 



^Pi.^.Pa^i, r 2 , r 3 + in) - W Pi ,p 2j p 3 (ti,T2,t 3 - in) 



W(t 3 + iir) — W(t 3 — in) 



|sin(27r6Pi)sin(^-2f|P 2 |)| Wi {Pi{ (t 3 , t 2 )Wi |p 2 |(r 3 , n) 



(4.65) 



b) when Pi E [-00, - |] and P 2 G [§, 00]: 



^1,^^3(71,72,73 + nr) - Wp u p 2i p 3 {fi,f 2 ,f 3 - in) 



W(f 3 + in) - W(f 3 - in) W Pi+ b P2 _t P3 (n,f 2 ,f 3 ) 



The cases bl), b2) and b3) are related to al), a2) and a3) by duality. 



(4.66) 



5 Microscopic formulation of 2D quantum gravity 
5.1 SOS model on a dynamically triangulated disc 



The functional integral of 2D Quantum gravity l|2.5j) is discretized by considering the statistical 
ensemble of all triangulations of the world manifold M . Simultaneously, the matter field defined 
on A4 should be discretized as a statistical system defined on such a dynamical triangulation |2j. 

Two-dimensional conformal field theories with central charge less than one can be constructed 
as the continuum limits of solvable two dimensional statistical models on regular lattices. There 
are several microscopic formulations of the c < 1 conformal field theories: as vertex models |4Uj . 
as height RSOS 41 and ADE models |42j . as loop models or as the (/-state Potts model or 
0(n) models |2()j . Correspondingly, we have several solvable microscopic formulations of the 2D 
quantum gravity with the matter central charge less than one. The statistical models are in this 
case formulated on a random lattice, the sum over geometries being replaced by the sum over 
all the lattice realizations. As a result we obtain solvable models of quantum gravity, which can 
be formulated also as matrix models |43 | ll5 [ l4"l] . 

All these models are mapped onto a special solid-on-solid (SOS) model with discrete un- 
bounded height variable h G 5Z and a Hamiltonian depending only on height difference of 
nearby sites. The map is called Coulomb gas picture [13110]. The SOS model in question de- 
pends on a continuous parameter g, the Coulomb gas coupling constant, and can be considered 
as a solvable microscopic realization of the gaussian field with background charge eo ~ g — 1. This 
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correspondence is very convenient since it allows to give the operator identification of the fields 
already of the microscopic level. Here we will restrict ourselves to the microscopic formulation 
of 2D QG based on the SOS model. 

We start with the definition of the SOS model on an arbitrary triangulated surface. (The 
case of lattice paved by squares has been considered in ^Ej.) The target space of the SOS model 
is an infinite chain, the Dynkin graph of Aao, whose nodes are labeled by half- integer heights 
h £ ~Z. The structure of the target space is determined by its adjacency matrix, A^, whose 
matrix elements are one if h and h' are nearest neighbors and zero otherwise: 

A hh> = 5 hh , + i+5 hh ,_i. (5.1) 

The adjacency matrix is diagonalized in the momentum space by the plane waves = e 2lnph 

i 

A hh , = [dp e 2 ™ v{h ~ hl) A p , Ap = 2cos7rp. (5.2) 
-l 

We will consider a non-unitary version of the SOS model, which is characterized by a background 
target-space momentum pq. Consider a triangulation T of the disc that is dual to a -trivalent 
planar graph. To each node i £ T we associate a height hi £ Z. The partition function of the 
SOS model on T is defined as the sum over all height configurations {i — * hi\i £ T}. 

On the flat lattice the Boltzmann weight of a height configuration is a product of factors 
associated with the the triangles 

( Sh \ 1/6 

Wa(/ii, h 2 , h 3 ) = T 5 hl h 2 5 h2h:i 5 hahl + 5 hlh2 A h2h:i A hahl J + cyclic, (5.3) 

where 

S h = e ~ 2 ™P° h (5.4) 

is the plane wave associated with the background momentum po £ [0,2]. The coupling T is 
usually called temperature. In presence of curvature defects one should also associate Boltzmann 
weights with the nodes j £ T: 

W.(hi) = (S h .)^/^ (5.5) 

where Rj is the curvature at the point j £ T. As we will see below, this factor is necessary to 
preserve the loop gas representation of the partition and correlation functions. 
The local curvature is defined as the deficit angle 

2-7T 

Ri = —(6- Ci ). (5.6) 

where Cj denotes the coordination number at the point i (the number of triangles having this 
point as a vertex). Similarly, the boundary curvature at the point j £ dT is defined as 

^ = |(3-9). (5.7) 

By the Euler relation the total curvature of the disc (including the curvature along the boundary) 
is 

Ri + 2 & = 47T - ( 5 - 8 ) 
ier iedT 
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The partition function of the SOS model on the triangulation T is given by the sum over all 
height configurations: 



Z ^= E U^MjlWA^h^hk). (5.9) 

{hi\i€T} i Aijk 

To define completely the partition function we should fix the boundary conditions. There are 
two natural choices of the boundary conditions in the height model. The fixed, or Dirichlet type 
boundary condition prescribes a constant height h along the boundary. The fixed boundary 
conditions where the hight takes a value a on the Dynkin graph are described by the type (a, 1) 
conformal boundary conditions [HS] 19 The free, or Neumann type boundary condition prescribes 
to sum over all heights along the boundary with weight one. 

The loop gas representation is obtained by expanding the partition function (|5.9j) as a sum 
of monomials, each monomial corresponding to choosing one of the four terms in (|5.3[) for each 
triangle. A neat geometrical interpretation of the monomials can be achieved if we represent 
graphically the r.h.s. of eqn. (|5.3|) as 




h, h, h, h, h, h, \ h 3 (5.10) 

In the last three terms the lines separate the two vertices with the same height from the vertex 
with different height. Then a monomial in the expansion of the partition function corresponds 
to a collection of closed non intersecting loops on the dual trivalent graph 7": 




The loops can be considered as domain walls separating the domains of constant height. The 
Dirichlet boundary condition means that the height of the domain adjacent to the boundary of 
T is set to the value hs- 

Let us consider the case of a Dirichlet boundary condition h = hs at the boundary. The 
weight of a loop configuration is obtained by summing over all possible domain heights. The 
Boltzmann weight of a domain T> with height h = h-p is the product of the Boltzmann factors 
associated with its nodes and triangles 20 . By the Euler relation (|5.8|) applied to the domain 
T>, the total power of Sh depends only on the number of boundaries of T>, which we denote by 
Nb(T>). As a consequence, the partition function ()5.9|) with Dirichlet boundary condition Kb 
can be reformulated as a sum over the domain configurations 

Z hB {T) = T^ £ T-*» £ U(s hv y~ NBiV) (5.11) 

loop configu- h-neiZ V 

rations on l £ 

19 In a regular lattice the general (a, r) boundary condition is realized by attaching r times fused weight to 
height a. This amounts to restricting the height variable to the interval [a — r/r, a + r/2]. The condition type 
(a, r) interpolates between the fixed (D) and free (N) boundary conditions. 

20 We will assume that a triangle belongs to given domain of constant height if at least two of its three vertices 
are inside the domain. 
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where £ tot is the total length of the loops, equal to the number of decorated triangles, and *4tot 
is the total number of triangles of T . We remind that in the sum the heights of neighboring 
domains differ by ±|. 

For any given domain configuration the sum over the heights can be performed using the 
fact that S x is an eigenvector of the adjacency matrix with eigenvalue n = 2cos7rpo 

^ A hh >S h i = 2 cos irp Q S h . (5.12) 
h> 

The summation starts with the innermost (simply connected) domains and it generates a factor 
n for each domain boundary. After the summation over all the domain heights, one obtains 
n* loops 5/ ls , where hs is the (fixed) height at the boundary. Thus we obtain a representation of 
the partition function as that of a gas of mutually- and self-avoiding loops on the triangulation 
T, with fugacity n = 2cospo an d tension logT 

Z hB (T) = S hB T Atot Yl T ~ hot (2cos7r Po )# lo °P s = S hg Z(T). (5.13) 

loop configu- 
rations on T 

If we fix, instead of the height, the momentum flowing through the boundary, we have 21 

Z ( P ) (T) =J2e- 2i ^Z hB (T) = S(P~Po) Z{T) 

h 

Thus the ground state of the SOS model is translationally invariant in the target space only if 
Po = 0; otherwise there is a non-zero momentum flowing through the boundary. 

The loop gas representation for Neumann boundary condition has been considered in |24j . 
The free boundary condition means that the height can change by ±^ at each point of the 
boundary. To define a height configuration we need to specify all heights along the boundary. 
In this case, the domain boundaries form not only closed loops, but also open oriented lines with 
ends on the boundary. There is an open line starting or ending at each boundary link. After 
summing over the heights using eqn. (|5.12j) . one obtains the partition function as the sum of all 
configurations of self- and mutually avoiding loops and lines. The closed loops acquire as before 
factors n = 2cos7rpo while the factors associated with the open lines ending at the boundary 
are 2 cos(7rpo/2). We can also consider Neumann and Dirichlet boundary conditions on different 
parts of the boundary. 

The microscopic theory of gravity is constructed as the statistical ensemble of all trian- 
gulations with given topology 2 . The trivalent planar graphs dual to triangulations can be 
interpreted as the Feynman graphs for the 't Hooft limit of a matrix model. In the case of the 
SOS and related models, such matrix model has been constructed in |44j . 

Consider first the ensemble of all triangulated discs 7^ made of A triangles and whose 
boundary has i edges. The disc SOS partition function with fixed height h at the boundary is 
given by the average of the matter partition function (|5.9|) 

W h (£,A) = M?) (5-14) 

21 A boundary state characterized by a fixed height is also referred as a Cardy state, while the state with fixed 
momentum p flowing through the boundary is called Ishibashi state. The two kinds of states are related by Fourier 
transformation |p> = J^he J- z e 2lriph \h). 
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where we assumed that the triangulations have a marked point at the boundary, so that there 
is no symmetry factor. Introducing the bare bulk and boundary cosmological constants A and 
z, we define the disc partition function with Dirichlet boundary condition h as 



<MA,,) = £ £ e-"-«™!&4. (5.15) 



Using the loop representation (|5.13|) . we can relate the disc amplitude with fixed height to the 
partition function of the loop gas on a random graph 

$ h (A,z) = e H *"> h $(A,z), (5.16) 



$(A,z) = Y.\ e- KA ~ zl Yl (2cosvrp ) [#loopsl e~ 2Me ^ (5.17) 

T loops on T 

where A is the area of the triangulation T, i is the length of its boundary and M = \ log T. 
The factor l/l takes into account the cyclic symmetry. 

5.2 Continuum limit and mapping to Liouville gravity 

In the case when T is the regular triangular lattice, the loop model Q5.1MJI has been solved by 
Bethe Ansatz technique |2Ul l40j . The loop model is critical in the interval < T < T c where its 
continuum limit is described by a conformal theory with central charge (we assume that po > 0) 

c dcnsc phaso = 1 -6pg/(l -po). (5-18) 

This is the dense or low-temperature phase of the loop gas. In this phase the loops occupy 
almost all the space, £ tot ~ «4tot- At the critical temperature 22 the loop gas exhibits a different 
critical behavior (the dilute phase) described by a larger central charge 

c d u utephase = l-6pg/(l +Po)- (5-19) 

In the dilute phase both the total length of the loops and the empty space diverge correspondingly 
as *4 e mpty ~ ^tot and £ tot ~ y/Atat- Above the critical temperature the loop gas is not critical, 
i.e. there is a finite correlation length determined by the typical size of the loops. A more detailed 
identification of the continuum limit has been made in |45|I20| . where it has been argued that 
the renormalization group trajectories of the SOS model flow to a gaussian field with electric 
charge at infinity. 

With Boltzmann weights defined by (|5.3jl - (|5.5|) this statement can be generalized to the case 
of a lattice with curvature defects ^H]- More precisely, the SOS model renormalized at large 
distances to the gaussian field (|2.2|) . with the identification 

e oX = 2np h. (5.20) 

Indeed, the effect of the curvature is to associate an additional weight factor 

exp(-^E J Rjvr/i i p ) (5.21) 

j 



The exact value of the critical temperature is T c = 2cos(7rpo/4) 
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with each height configuration. The logarithm of this factor represents a discretization of the 
linear term of the gaussian action Q2.2|) . which leads to the identification (|5.2U|) . Prom (|5.18|) 
and 1)5. 19J) one gets, assuming that < b < 1, 

e o = Po/b, X = 27r6/i in the dense phase; 

eo = bpo, x = 27i7i/£> in the dilute phase. (5.22) 

Thus the gaussian field x with background charge en. gives the effective field theory for the 
large-distance behavior of the SOS model with background target-space momentum po, with the 
identifications (|5~2TIj) and (|5~2^|) . 

In the same way as all c < 1 conformal field theories are mapped onto the gaussian field, one 
can map the underlying solvable spin models (ADE, 0(n), Potts, ...) onto the SOS model in the 
sense that the partition and correlation functions can be interpreted in terms of distributions 
of electric and magnetic charges. This mapping is known as Coulomb gas picture |45[l20j. 
Traditionally one uses the Coulomb gas coupling constant 

g = 1 — po = b 2 in the dense phase; 

g = 1 + po = l/b 2 in the dilute phase. (5.23) 

Then (|5.22|) reads in both phases 

eo = Po/Vg, X = 2vr^5 h. (5.24) 

The SOS model and related models on a random graph have been solved in JH] (for the 
dense phase) and in ^H] (for the dilute and the multicritical phases). The bulk critical behavior 
due to triangulations with divergent area is achieved along a line in the two-dimensional space 
of coupling constants {A, M} ending at a tricritical point {A*,M*}. The critical line describes 
the dense phase of the SOS model coupled to gravity, and the tricritical point - the dilute phase. 

From now on we will consider only the scaling regime and will denote by the same symbols A, 
M and z correspondingly the renormalized cosmological constant, loop tension and boundary 
cosmological constant, correspondingly. The scaling of the bulk and boundary cosmological 
constants with the renormalized loop tension 2M is |18U24j 



A dilute ^ M 2^ A dense ^ M 2g 

^Diriehlet _ M) z Neumann_ M <; (in both phases). (5.25) 

The Dirichlet boundary can be considered itself as a domain wall, which explains why ^Dmchlet 
scales as the domain wall tension. The scaling (|5.25j) leads to the following identification of the 
bulk and boundary Liouville interactions: 



Bulk Dirichlet boundary Neumann boundary 

Dilute phase: g = l/b 2 > 1 : fx e 2 ^ fi B & H e& b ^' 2 ^ 

Dense phase : g = b 2 < 1 : \i e 2 ^ fl B \i B e b<t> 
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Each value of the central charge in the interval —2 < c < 1 is realized both in the dense and 
in the dilute phase of the SOS model. The phase that connects to the continuum path integral 
is the dilute phase with Dirichlet boundary or the dense phase with Neumann boundary. The 
dilute and dense phases are related by a T-duality transformation for the matter field. In order 
to have identical conventions with the continuous approach, we normalize the bulk cosmological 
constant so that A dilute = M 2 and A dense = M 2g . 

The correlation functions of the boundary operators are meromorphic with respect to the 
boundary cosmological constants z, with a cut [—00, — M] along the real axis. When the bound- 
ary cosmological constant approaches the value — M, the partition function diverges because of 
the dominance of world sheets with infinite boundary length or, in terms of the Liouville path 
integral, because of the diverging integral with respect to the zero mode of the Liouville field at 
the boundary. The branch point is resolved by by an uniformization map 

z = Mcos\it. (5.27) 

The complex z-plane is mapped to the strip |Imr| < n with the points r and — r identified. The 
two sides of the strip parametrize the two edges of the cut. The branch point corresponds to 
the points r = m = —in. 

5.3 Bulk and boundary operators in the SOS model 



Bulk electric and magnetic operators: The electric charge in the SOS model corresponds 
to the target space momentum p. Since the momentum space is periodic, p + 2 = p, we can 
restrict p £ [0, 2]. The electric charge is carried by the electric, or vertex, field 

V {p \j) = e ^iPo-p)hj. (5.28) 

In terms of the loop gas the electric operator (|5.28fl changes the weight of the non-contractible 
loops, i.e. those encircling the point where the operator is inserted. Instead of n = 2cos7rpo, 
each non-contractible loop gets a factor 2cos-7rp 20 . It follows from the matching of charges 
(|5.22|) that the bulk dimensions of the SOS vertex operators (|5.28f) are 

A y(rt = t^A. (5.29) 

Unlike the electric operators, the magnetic, or vortex, operators S^ L ' have infinite but discrete 
spectrum L £ Z. In terms of the loop gas these operators are the "star polymers" |21[ I46|. 
representing configurations of L open lines starting at the same point. Strictly speaking, the 
magnetic operator S^ L ' on the lattice is associated with a hole with boundary of length L, with 
an open line starting at each edge. For example, the configurations associated with the operator 
look as follows: 
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Since in the SOS model each line represents a domain wall separating two neighboring heights, 
the magnetic operator creates a a discontinuity 5h around the small loop, which can take values 
Sh = —jg,—72 + l,...,7£. When 5h = ±-|, the operator can be identified with the spinless 
magnetic operator O e0tm for the gaussian field, with m = ±L^/g/4 and conformal dimension 23 



vl 



±9 



Al/2.o in the dilute phase (g > 1) 
Ao,L/2 in the dense phase (g < 1) 



(5.30) 



Boundary electric operators: The boundary vertex operator is defined by inserting the 
vertex field 



V P (j) 



o 2-Ki(\p -p)hj 



(5.31) 



at a point on the boundary with Neumann boundary conditions on both sides. In terms of the 
loop gas the boundary vertex operator changes the weight of the noncontractible lines, i.e. the 
lines that connect the boundary on its left and the right. These loops get additional weight 
factor cos(7rp)/ cos(7rpo/2). Its conformal dimension follows from 1)2.16)) : 



V - vj 

4<7- 



(5.32) 



Boundary magnetic operators: The boundary star operator Sl creates L lines at some 
point of the boundary with Dirichlet boundary conditions: 

L lines 




The operator Sl implies a discontinuity of the boundary height 5h = bright — ^ieft> which can 
take values 5h = + l,---, Again, when 5h = ±-§ the operator Sl can be identified 

with the boundary magnetic operator for the gaussian field with magnetic charge m = ±LyJ~g/2. 
Therefore its conformal dimension is given by (j2.17[) : 



As, 



(Lg - p ) 



Pi 



±9 



Al+i,i in the dilute phase (g > 1) 
Ai 5 i + i in the dense phase (g < 1). 



(5.33) 



23 In the rational points b 2 — p/(p + 1) the star operators correspond to the magnetic operators 
the minimal models, with JV = . 



>N-l-n,N 



of 
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Boundary twist operator: The boundary twist operators are associated with the points 
separating the Dirichlet and Neumann type boundaries: 




To our knowledge, the boundary twist operator operator has not been discussed in the 
context of the loop gas on a flat lattice. It seems to have a simple description in terms of the 
gaussian field only for 6 = 3/4 and 6=1. For generic 6 we found its flat scaling dimension from 
its gravitational dimension using the KPZ scaling formula: 



A 



T 



(3 - 26 2 )(26 2 -l) 
166 2 



(5.34) 



Thus the boundary twist operator is not a degenerate field, but can be classified as a "half- 



degenerate" field with r 



2' 



s = 0. Note however that in the case of a minimal model 6 = — tt 



? +1 

with p odd, the boundary twist operator can be identified with the diagonal operator {r, s\ = 
{p + l,p + 1} of the Kac table. For example, in the Ising model (6 = ^3/4) this is the order 
operator a (r = 2, s = 2) with dimension 1/16. 

We will also consider excited boundary twist operators Tl obtained by the fusion of a twist 
operator T and a star operator Sl- 




The dimension of such operator is 



(L + l) 2 /6 2 -(6-l/6) 2 
t l = j = A L+ i . (5.35) 



Here again we extracted the conformal dimensions from the two-point correlation function in 
2D QG using the KPZ relation. 
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6 Difference equations for the boundary correlators in the mi- 
croscopic approach 

6.1 Factorization property of the sum over geometries 

The boundary correlators in the SOS model satisfy a set of loop equations, which can be restated 
as functional equations for the boundary r-parameters. These equations are consequence of the 
remarkable factorization property of the functional measure of the loop gas on a random surface. 

Let r be a self-avoiding line on the disc M connecting two points on the boundary. It splits 
the disc M. into two pieces M.\ and M 2 such that M. = M.\ U M. 2 . We can first perform the 
integral with respect to all configurations of the self- avoiding line T, and then integrate with 
respect to all metrics g a b on A4. Alternatively, cutting open the path integral along the line 
r, we can integrate first with respect to the metrics on A4± and M 2 and then with respect to 
the configurations of the separating curve T. Due to the general covariance of the measure, the 
result of the integration depends only on the length £ of the curve T. Therefore the path integral 
over the configurations of the curve T reduces to one-dimensional integral along its length i. 
Thus we have the following factorization formula 

J Vg ab (M) jw = P ' d£Vg ah (M x ) j Vg ab (M 2 ) (6.1) 
M ° 

where DT is the measure over self-avoiding lines r E M with fixed endpoints. This formula 
is quite obvious from the microscopic realization of the measure as a sum over triangulations 
(see US] for details). 

Now we use the fact that the SOS partition function on A4 in presence of such line factorizes 
into a product of two SOS partition functions on .Mi and on M 2 with Dirichlet boundary 
conditions along C: 

Z{Mi UM 2 ) =Z(Mi)Z(M 2 ). (6.2) 

As a consequence, we can represent the partition function of the SOS model in presence of a 
domain wall as an integral with respect to the length of the domain of the product of the SOS 
partition functions, for which the domain has the meaning of a Dirichlet boundary. 

6.2 Disc loop amplitude 

The most basic observable is the disc loop amplitude Wh (£) , defined as the partition function of 
the loop gas on a random surface with the topology of a disc with fixed height h on the boundary 
(Dirichlet boundary condition), boundary length I and a marked point at the boundary 
Since we are concerned only by the scaling limit, by length we understand the renormalized 
length of loops or Dirichlet boundaries. The length is related to the constant mode of the 
Liouville field on the boundary as 

| e b ^° in the dilute phase, 
[ e^l h in the dense phase. 

It follows from eq. Q5.11JI that W} t {€) and its Laplace transform 

W(z) =f^°d£e- ze W(e) = ^^ (6.3) 
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depend on the boundary height as 



W h (i) = S h W(£), W h (z) = S h W{z), 

(where W{£) and W{z) are the loop gas disc partition functions with a marked point at the 
boundary). 

The disc amplitude satisfies a quadratic equation of motion [T61IT8] . Here we give a heuristic 
derivation of this equation in the continuum limit. Consider first the ensemble of empty world 
sheets, called pure gravity. It corresponds to the values g = 1/2 (c = —2) and g = 3/2 (c = 0) 
of the SOS model. It is known that in this case the loop amplitude W(z) satisfies a quadratic 
equation 
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W 2 (z) 



where P(z) is a quadratic polynomial in z, whose coefficients depend on the bare couplings Mq 
and fj,Q. Written in the ^-representation 25 



d£' W(£')W{£ - £' 







(*>0) 



(6.4) 



this equation has the following meaning: a small variation of the world sheet metric at the 
marked point produces a contact term associated with the degenerate world-sheet geometries, 
for which another point of the boundary is at zero distance from the marked point. The integral 
can be regarded as the total contribution of these degenerate geometries, for which the integral 
over metrics factorizes into a product of two disc amplitudes. We can write (|6.4jl symbolically 
as 








(6.5) 



Now consider the general case of matter coupled to gravity. Apart from the contact term 
that comes from the degenerate geometries, there is a second contact term associated with a 
loop coming close to the marked point 26 . By the factorization property we mentioned above, 
the contribution of such loops is given by an integral over their length times the loop fugacity 
n = —2 cos 7rg 



f 

Jo 



d£'W{£')W{£ - £' 



or pictorially 



2cos7T# / d£' W(£')W (£ + £') 
/o 



= 







(£ > 



(6.6) 





(6.7) 

The Laplace-transformed loop amplitude VF(2:) is a rational function of z with a cut extending 
between the simple branch point at z = —M and a branch point of (generically) infinite order 
at infinity. The position of the right branch point is related to the tension at the boundary, 



W{£) 



-Ml 



It is consistent to assume that the tension at the boundary is equal to half of 



24 We have shifted W(z) by a polynomial, which is irrelevant in the continuum limit. 

25 Here we require that £ is strictly positive, because we are not interested in the singular at I = terms. 

26 This contact term is in fact the boundary one-point function of the energy operator of the loop gas. The 
energy operator is primary in the bulk, but appears as the descendent of the identity operator when classified 
according to the irreducible representations of the boundary Virasoro algebra |29|. 
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the tension of the loops in the bulk. Then weight factor e~ 2Mi of the loop in the second term 
of 1)6. 5|) gives the correct boundary tension for the loop amplitudes associated with the internal 
and the external domains. In terms of boundary cosmological constants, eq. 1)6.5)1 reads 27 



W\z) — 2 cos ng / 

2m 



dz' W(z)-W(z') 



z — z' 



W(-z') = 



(6.8) 



where the contour of integration encircles the cut — oo < z < —M. Then equation 1)6. 8|) implies 
the following condition relating the real and imaginary parts of VF(^) along the real axis 



ImW(z) [W(z + iO) + W(z — iO) — 2 cos ng W{-z)] = 0, (z £ 



(6.9) 



In the parametrization z = Mcoshr, which unfolds the branch point at z = —M, the 
upper and the lower sides of the cut correspond to the lines r + in and r — in, < r < oo. 
Therefore 16.91 can be written as the following finite-difference equation for the analytic function 

W(t) = W(z{t)) 

W(t + in) + W(t -in) -2 cos ng W(t) = (6.10) 
which is solved up to a constant factor by 



W(t) 



1 



sin ng 



M 9 cosh(gr). 



(6.11) 



6.3 Bulk one-point function 

The bulk one-point function (z) is the partition function of the loop gas on the punctured 
disc with a vertex operator V p inserted at the puncture and Dirichlet boundary condition char- 
acterized by a height h and boundary cosmological constant z = Mcoshr. It has the form 



U^(z) 



where (z) is the loop gas partition function on the disc in which the fugacity of the loops 
encircling the puncture is changed to 2cos7rp )16j . In order to write the loop equation, we 
should first differentiate with respect to z in order to create a marked point on the boundary. 
The derivative 

R(P)( Z ) =dU M {z)/8z, (6.12) 
which is also a bulk-boundary correlation function, satisfies the following integral equation ^B] 

2R(r\z)W(z) + 2COS7TPQ / — — ^— ^V (-z') 

J 2ni z — z 

J 2m z — z 

or pictorially 




+ 




W t R * + 




+ 








(6-14) 

27 The integral is divergent at infinity, but we assume that the amplitude W(z) is regularized so that the cut 
ends at a second branch point at large but finite distance. The integration then can be done using the Cauchy 
formula, and then the cutoff may be lifted. The rigorous procedure is explained in |16| and |18|. 
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From here we obtain, for z G 



lmW{z)[R^ p \z + iO) + R (p) {z - *0) +2 cos 7rp R^ p \— z)\ 
+IolRW + iO) + VF(z - iO) + 2 cos vrp = 0. 



(6.15) 



The second line of vanishes due to Q6.1UJI and one finally obtains the following functional equation 
in the r-parametrization 



(r + in) + R {p) (r -in) + 2 cos np flW (r) = 0. 



The solution of this equation with the correct scaling is given by 

R(P)( T ) =d z U® ~ in- 



_| p |_ 1 sinhpT 
sinhr 



(6.16) 



(6.17) 



This expression is to be compared with the expression (|3,27|) found from Liouville theory. Inte- 
grating with respect to z, we find 



U {p \t) ~ Ml p l coshpr. 



(6.18) 



in accordance with Q3.14JI . 

One can easily write loop equations like (|fi.5|) and (|6.14j) for any boundary correlation function 
of electric and magnetic operators under the condition that there should be at least one with 
Dirichlet boundary condition for the SOS field. The r.h.s. will contain a contact term due to 
a pinching of the disc, another term associated with the energy operator at the boundary, and 
a number of contact terms, each associated with a line or a non-contractible loop touching the 
marked point. 

In the next section we will consider another class of functional equations, which follow only 
from the factorization property of the measure. These equations will be shown to be equivalent 
to the difference equations in Liouville gravity. 

6.4 The boundary two-point function of magnetic operators 

In this subsection we recall the calculation of the boundary two-point function of magnetic, 
or star, operators presented in j2^j. To begin with, consider the simplest example of the disc 
amplitude with two marked points and two different boundary cosmological constants, 



W( Zl ,z 2 ) 



di e 



-Z\l\—Z2l2 



Wih+h) 



W( Zl ) - W(z 2 ) 



z\ - z 2 



(6.19) 



This is a particular case (L = 0) of the correlation function of two gravitationally dressed star 
operators 

D L (Ti,r 2 ) = (S T L ^SZ 1 ). (6-20) 

In the loop gas representation, this correlation function is the loop gas partition function in 
presence of L open non-intersecting lines connecting two points on the boundary: 




- X, 




(6.21) 
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The inverse Laplace transform of Dl, which is the partition function for fixed lengths £\ and £2 
of the two boundaries 



D L (n,T 2 ) = / dhd£ 2 e 
lo 



-z 1 e 1 -z 2 £2 



D L (£ 



1,^2 



(6.22) 



can be expressed through the disc partition function W(£) as follows. Consider the sum over 
all loop configurations with fixed lengths £[, ...,£' L of the branches of the star polymers, which 
is given by the product of L + 1 disc partition functions, and then integrate with respect to the 
lengths 



n L {h,£ 2 ) 



W(h + £' X )W{£\ + £' 2 )...W{£' L + £2). 



(6.23) 



The integral representation (|6.23j) is equivalent to the recurrence relation between Dl and 
Dl-i 





D L (h,£ 2 )= d£W{h+l)D L ^{£,£ 2 ), 
'0 



D, 



0^1,^2) 



W(£ x + £ 2 ) 



or, after a Laplace transformation, 

D L (z 1 ,z 2 ) = 



dz W{zi)-W{z) 



2ttz 



Z\ — z 



D L -i(-z,z 2 ) 



(6.24) 
(6.25) 

(6.26) 



where the contour of integration encircles the cut —00 < z < —M. This integral equation 
implies a condition on the discontinuity of Dl(z) along the cut: 

lmD L (z,z 2 ) =ImW(z)-D L _ 1 (-z,z 2 ), z < -M. (6.27) 

In terms of the variable r, this condition is equivalent to a finite-difference equation |25| 

D L (n +m,T 2 ) -D L (n -iir,T 2 ) = [W{ti + m) -W{n - m)]D L - l {T l ,T 2 ) (6.28) 

or, taking into account the explicit expression (|6.11j) for W(t) 

smiird/dn) D l (t 1 ,t 2 ) = M 9 sinhgn D l _i(t 1: t 2 ). (6.29) 



Together with the initial condition Dq(t\,T2) = W(t\,t 2 ), eq. (|6.19j) . this recurrence equation 
has the unique solution 



D L (zi,z 2 ) 



CL 



Zl ~ {~) L Z2 



\\ (W(n + iir(L - 2k)) - W{t 2 )) 



(6.30) 



fc=0 



where the normalization constants are given by 

L 



CL 



h 1 n 



sin irg 



^ sin(fc + l)irg 
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It is a rational function of the variables z a = Mcoshr a and W a = M 9 cosh gr a (a = 1,2): 
Wi -W 2 



D (zi,z 2 ) 
D 1 (z 1 ,z 2 ) 

D 2m {zi,z 2 ) 

D 2m+1 (zi,z 2 ) 



Z\ - z 2 
ci 
Z\ + z 2 

W- 

c 2m 



Wl + Wf - 2 cos itg WiW 2 - M 29 ) 



c 2m +i 
z\ + z 2 



^» T](w? + W 2 - 2cos2^ Wl W 2 - M 29 ^^-\ 
i - z z £Ji V ' sin 7T# y 

IT (V + Wi - 2cos(2A: + 1)tt 5 W 1 W 2 - M 2 ° S ^ 2 \ + 1)tt<A 

ii V _ sin 7T ff y 



Let us compare (|6.28[1 with the difference equation obtained in Liouville gravity. First, notice 
that this requires the choice (|4,6U|) of the constant k. 

a) g > 1. 

Then the operator Sl should be identified with the boundary field £>l+i,i with Pl+i,i = 
^eo + L^g. Then \P\ > ^ and eq. (j6.28j) reproduces the difference equation (|4,62[) in Liouville 
gravity 

b) \ <g<l. 

Then the operator Sl should be identified with the boundary field £>i,l+i with P^l+i = 
|eo — L|. Then |P| > | for all L > 1 and eq. (|6.28f) coincides with the difference equation 

In the interval < g < ^ the statistical interpretation of the loop gas partition function does 
not exists because the fugacity of the loops is negative: — 2cos(7r<7) < 0. 

6.5 Boundary three-point functions of magnetic operators 

The same techniques can be used to write difference equations for correlators of more than two 
operators. Consider the correlation function of three boundary magnetic operators Sl 1 , Sl 2 
and Sl 3 

C Ll ,L 2 ,L 3 (n,r 2 ,r 3 ) = (S^S^S-f ) (6.31) 

or pictorially 



_ x 





(6.32) 

Let us first evaluate the three-point function of the identity boundary operator with L = 0. 
This is the partition function on the disc with three different cosmological constants Z\ , z 2 , z% 
along three segments on the boundary. It is equal to the integral 



C 000 (z 1 ,z 2 ,z 3 )= / dt x ...dh e-^^Wih+h + h) 
Jo 



dz W(z) 
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or, after performing the integration, 

n ( v (*2 ~ Z 3 )Wi + (Z3 - z x )W 2 + (zi - Z 2 )W 3 ^„„v 

^000 Ulj %2, Z 3 ) = -. r~. -7 r . (0.66) 

(Zl - z 2 )(z 2 - Z 3 )(z 3 - Zl) 

This expression solves the Liouville difference equation (|3.25|) and coincides, up to a multiplica- 
tive constant, with the Liouville three-point function ([3.28)1 . 

It is straightforward to write a recurrence equation relating the correlation functions Cl 1 l 2 Lz 
and Cl 1 _i.l 2 _i.l 3 by using the factorization property associated with the most external line 
relating the points 1 and 2: 




L 2-! (6.34) 

Skipping the derivation, which is the same as in the case of the two-point function, we write 

Cl 1i l 2i l 3 (t3 + iir,Ti,T 2 ) - Cl u l 2 ,l 3 {tz - wt,tl,t 2 ) = 
[W(t 3 + in) - W(t 3 - in)} C Li ^ 1m . 1m (t 3 , n , t 2 ). (6.35) 

Knowing the explicit form of the function for Cooo > given by eq. 1)6.33)1 , we can give the following 
Ansatz for the general form of the solution 

r T. k z k P k LlML3 (w^w 2 ,w 3 ) 

Clm " = rw* -(-i)—-,) (6 - 36) 

where L\ = n 2 + n 3 ,L 2 = n 3 + ni,L 3 = n\ + n 2 and Pli,l 2 ,l 3 is a polynomial of degree 
i(Li + L 2 + L 3 ) + 1 = m + n 2 + n 3 + 1. 

Let us compare eqn. ()6.35)1 with the difference equations for the three-point function obtained 
in Liouville gravity in sect. 4: 

a) g > 1. 

Then the operator Sl should be identified with the boundary field Bl+i,i with Pl+i,i = + 
L^g. In this case eq. ()6.35)) coincides with the difference equation 1)4.62)1 and \P\ > ^. 

a) \ < g < 1. 

Then the operator Sl should be identified with the boundary field £>i,l+i with = 
In this case eq. ()6.28)) coincides with the difference equation 1)4.66)1 and \P\ > |. 

Remark 1: Eqn. ()6.35|) was derived in the SOS model, where the loop amplitude W(r) does 
not depend on the height of the boundary. For a general ADE model this equation is written 
as follows. Let C ' L l '_i l 2 -i l 3 -( Ti ' T2 ' r 3) ^ e * ne boundary 3pt function with Liouville and matter 
boundary conditions a k , h k (k = 1, 2, 3). Then the difference equation for (3 k = Q— [(Lk + l)g— 1] 
takes the form 

^;KL 3 3 (n,T 2 ,T 3 + l n) - ^;t,L 3 3 (n,r 2 ,r 3 - zvr) = 
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= [W h3 (r 3 + «r) - W h3 (r 3 - «r)] A hsh , <?£:5;S-i,Xsfa.^) (6.37) 

where W^(r) = S^M 9 cosh gr, 5^ is the Perron- Frobenius vector and A^/ is the adjacency 
matrix of the Dynkin graph. Writing the solution in a factorized form 

C^U,lM^) = ^-l 2 ;S-i >i3 C '^-^-i^(n,r 2 ,r3) 
we get eqn. (|6.35j) as well as another equation for the matter factor. 

Remark 2: We find it interesting to propose the following formal representation of the three- 
point function as the integrated product of three two-point functions and one loop amplitude 




(6.38) 

, 3 _i(-z,z 3 ). (6.39) 



where n± = 5(^1 + L2 — L3), 122 = 2(^1 + L3 — L2), = 2(^2 + L3 — L\). Written in terms 
of r- variables, 

M 9+1 f°° 

Cxi,l 2 ,l 3 (ti ) T2 ) T3) = — / drsinhrsinhstT L> ni _i(r,ri)Z) n2 _i(r, r 2 )L» n3 _i(r, r 3 ). (6.40) 

2ir J_ OQ 

This expression satisfies the difference equations on each of the variables ti,t 2 ,T3 and for L\ = 
L2 = L3 = (i.e. n\ = 77-2 = ri3 = —1) it reduces to the expression ()6.33l) . 

Although the three-point function in Liouville theory has an integral representation ()3.23j) . 
this seems very different from (|6,4()jl . In particular, we were not able to evaluate directly the 
three point function in LFT from its integral representation (|3.23|) : we managed to give 
its dependence with respect to the boundary parameters onlythanks to the functional equation 
(|3.25|) . We find it remarkable that its equivalent in 2D gravity is so simple to compute. 

It is natural to expect that eq. (|6.39|) has its counterpart in Liouville theory, namely the 
difference equation that follows from the fusion rules with a higher degenerate field B\^ + i or 

6.6 Relation between the bulk-boundary and the boundary-boundary two- 
point function of star operators 

Let Rl{z) be the structure constant for the fusion of a bulk and boundary L-star operators 
R L (r) = (S^ S l >: 




(6.41) 
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The factorization of the measure applied for the union of the two most external lines leads to 
the identity 





which reads, in ^-representation, 

R L (z) = jdz'W{z,z') D L ^ 2 (-z',-z'). 

The condition on the cut 

ImR(z) =ImW(z) D L - 2 (-z, -z), z < -M 
gives the following finite-difference equation 

R l (t + in) - R l (t - m) = [W(t + iir) - W(t - m)) D l ^ 2 (t, t) 



(6.42) 



(6.43) 



(6.44) 



The function Rl(t) should be proportional to the Liouville bulk-boundary structure constant 
R a {a,(3) with a = Q/2-L/4b and = Q/2- (L + l)/2b + b/2 = b-L/2b. Indeed, with the help 
of the identity IJH.22J) . eq. (|6.44j) can be identified as a particular case of the Liouville functional 
equation (|3.2U|) . 



6.7 The boundary two-point function of twist operators 

We have checked that the functional equations in the SOS string and in Liouville theory are 
compatible for two infinite discrete sets of equally spaced Liouville charges. Now we are going to 
consider the case when some parts of the boundary have Neumann boundary condition. This will 
allow to compare the two approaches for situations involving a continuous spectrum of Liouville 
charges for some of the operators. 

The simplest correlation function involving both Dirichlet and Neumann boundary conditions 
has been calculated in |2l] . This is the correlation function Q(z, z) of two boundary changing 
operators that intertwine between the Dirichlet and Neumann boundary conditions (see Fig. 10 
of ref. |2l] ). Recently it has been checked [22] that this correlation function coincide (again up 
to a numerical factor) with the boundary Liouville correlation function ()3.16j) for 

P = b/2 + 1/46. (6.45) 

For completeness we review here the observation of |25j . We denote by 

z = McosIit, z = M 9 coshgf (6.46) 

the cosmological constants along the Dirichlet and Neumann type boundaries, correspondingly 




(6.47) 



42 



Applying the factorization property to the maximal chords we get the following integral equation 





(6.48) 

whose simplest analytic form is in the (£, z) representation, where I is the length of the Dirichlet 
boundary and z is the cosmological constant associated with the Neumann boundary: 

(l(l,z) = 1 + Vn n / d£t...d£ n W(£ + £ 1 + ... + £ n ,z) (l{l 1 ,z)..h(t n ,z% (6.49) 

n=0 J ° 



where n = 2cos(|vrpo) = 2sin(7rg/2) is the fugacity of the open lines. From here we find for the 
Laplace transform 

Q(z, z) = l- J °° d£e- ez n(£, z) 

1 f dz' 1 W(z') 

n(z,z) = ; —- & —, \ \ . (6.50) 

v ; 2sin(7T5/2) J 2iri z - z' n(-z',z) v ; 

Taking the imaginary part along the cut, we get the functional equation of [21] 

Imn(z,~z) = j—- ^ W(Z \ (6.51) 

v ' 1 2sin(7r 5 /2) n(-z,z) v ; 



or, in terms of r, 



0(T + ^7^,T)-0(r-^7^,T) = — — — . (6.52) 

2 sm(irg/2) ll(r, t) 



Let us compare eq. (|6.52jl with the difference equation 1)4.55)1 for the boundary two-point 
function in Liouville gravity. The two equations are compatible for P = i and 

n(r,f) = ^= Wi (r,f). 
V2 3^ 

This confirms the identification of the gravitationally dressed twist operator as the boundary 
state Bp with P = i: 

T = (z6) 1 / 2 2~ 1//4 £>j_. (6.53) 

4b 

This is a "half-degenerate" boundary field with r = i and s = and gravitational conformal 
dimension 

2b 2 - 1 

The flat dimension of this operator is given by 1)5.34 J) . 

We can consider, more generally, the excited boundary twist operators defined as a juxta- 
position of a boundary twist and magnetic operators 

T L = S L T. (6.54) 
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Geometrically the operator is constructed as the source of L nonintersecting lines inserted 
between at the point where Dirichlet and Neumann boundaries meet. The two-point function 

n L (r,f) = (TfTY) 

satisfies the same recurrence equation ()6.26j) 

Q l (t + in, f) -n L (r- in, f) = [W(t + in) - W(t - m)]J2i_i(r, f ), (6.55) 

but with different initial condition VLq(t,t) = 0,(t,t). Comparing with the difference equations 
in Liouville quantum gravity we identify the excited twist operators as fields "half-degenerate" 
boundary KPZ fields r = L + ±,s = with P = (L + ±) /b - b. 

6.8 The boundary correlation function of two twist operators and one vertex 
operator. 

Here we will consider the correlation function of two twist operators and an electric operator 
with arbitrary target space momentum p: 

%(t, h\h,f 2 ) = (T Tfl V p flf2 T f2r ). (6.56) 

In this case the correlation function depends explicitly on the height h of the Dirichlet piece of 
the boundary between the two twist operators. The corresponding pattern of domain lines and 
loops is: 




(6.57) 



There are two kinds of open lines: the contractible ones whose ends belong to the same Neumann 
boundary and the non-contractible ones that connect two different Neumann boundaries. The 
non-contractible lines have additional weight factor cos7rp/ cos-7r(^po)- 

Due to the momentum introduced by the vertex operator V p the correlation function has 
nontrivial dependence on the height h: 

V p (T,h\n,f 2 ) = e 2mph %(r\f 1 ,f 2 ). (6.58) 

The function ^? p (z, h\z\, z 2 ) satisfies a linear integral equation, obtained by factorizing with 
respect to the non-contractible line closest to the Dirichlet boundary. The equation involves 
the 4-point function T(z, z'\z\z 2 ), which is the partition function on a disc whose boundary is 
divided into four segments, with alternating Dirichlet and Neumann boundary conditions and 
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no open lines connecting the two Neumann boundaries: 28 
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(6.59) 



We can write 1)6.57(1 as the sum of a term containing no non-contractible lines, which we denote 
by VP 1/2, and a term that factorizes into T and fy p : 




+ 



r 









(6.60) 



This yields the integral equation 

—X{z,z'\ziZ2)^p{z'\ziZ2)- (6.61) 

The 3-point function vPi/^tIti, t 2 ) is defined in the same way as Vtp, but since cos(|-7r) = 
there are no lines connecting the sides f\ and r 2 . It is not difficult to see that the function T is 
expressed through ^>\i 2 as 



T(z, z'\z\Z 2 ) 



VP 1/2 (z|5i, z 2 ) - ^ 1 / 2 (z'\z 1 ,z 2 ) 



z — z' 



(6.62) 



The 3-point function ^i/ 2 itself is expressed in terms of the two-point correlator of twist operators 
fi(*,z) (see sect. 6 of EH ) 



d( W{z) - W(() 



[2cOs(i7Tp )] 2 J 27Ti Z-C fi(C,*lMC,*2)" 



(6.63) 



Again integral equations (|6.61j) and (|6.63|) can be turned into finite-difference equations for the 
r- variables: 



sin7r<9 T vP p (r|fi,f 2 ) = [sin7r<9 r ^i/ 2 (r\fi, f 2 )\ [1 + cos7rp VP p (r|fi, f 2 )] 



(6.64) 



sin7T<9 T * 1/2 (T|fi,f 2 ) 



sin gr 



[2cos(ivrpo)] 2 ^(T,n)0(r,f 2 )- 



(6.65) 



28 The function T is actually the correlation function of four boundary twist operators, T(z,z'\ziZ2) = 
^ T j*z2 T « 2 z' T *'2i T 2i*<} ; with target-space momenta p = 1/2 and p = -1/2 associated with the two Dirich- 
let boundaries (with boundary parameters r and r'). The effect of inserting the momentum 1/2 is that all lines 
connecting the two Neumann boundaries will acquire a factor 2cos(^7r) = 0. 
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The non-homogeneous term on the rhs can be eliminated by shifting ^ p —* ^ p — 1/ cos7rp. As 
a result we finally arrive at the homogeneous equation 



cos 7rp smgr 

smird T %(t\ti,t 2 ) = -rn — of _ - ; * P (r n,^) (6.66) 



This equation cannot be compared directly to the difference equations obtained in Liouville 
quantum gravity, which were written for too special situation (different signs of the two adja- 
cent momenta and the neutrality condition imposed.) Nevertheless is we consider the formal 
expression 

([B 1/4b r[B bp f f [B 1/4b f T ) disc , (6.67) 

the corresponding difference equation 1)4.65 J) is almost identical with ()6.66)) . The extra factor 
sin(7rg/2) = sin(7r/26 2 ) in the denominator can be explained with the fact that here both fields 
have positive momentum P = 1/26. The factor cos 717? has its origin in the matter field. In the 
second term on the r.h.s. of eq. 1)6.60)) the height of the Dirichlet boundary of is either h + \ 
or h — ^, because it is separated by one domain wall from the height h. The sum over the two 
possibilities gives, taking account of 1)6.58)1 . the factor 2cos7rp. 

Note that we can replace the twist operators T by excited twist operators Tl = S^T and 
apply the factorization condition to one of the lines of the star operators. This will produce eq. 
dQ2t with Pi = (L + i)^ - \ and P - 3 = bp for all L > 1. 



7 Conclusions 

In this paper showed that the the continuous and the discrete approaches to 2D QG lead to 
the same boundary correlation functions. First we considered pure Liouville theory and showed 
that all Liouville boundary structure constants satisfy linear difference equations with respect 
to the boundary parameters. 

We interpreted these equations in the context of Liouville quantum gravity with gaussian 
matter field. We observed that after rescaling the boundary fields by the leg factors and taking 
into account the charge neutrality condition, the charge dependence of the difference operator 
disappears. The difference equations have the form of recurrence relations for the charges. After 
repeated application of the difference operator, one obtains a correlation function in which some 
of the matter fields have "wrong" Liouville dressing. If this is the case, we get, using the 
reflection property, nonlinear difference equations for the physical fields. 

Then we considered the microscopic realization of 2D QG and defined the boundary corre- 
lators as expectation values in the ensemble of non-intersecting loops and lines on a randomly 
triangulated disc. We derived difference equations of the same form using the factorization of 
the measure over random surfaces after "cutting it open" along a non-intersecting line on the 
world sheet. 

Comparing the difference equations obtained in the two approaches we see that they become 
identical for certain normalization of the boundary one-point function. We consider this as 
an important "experimental" confirmation of the Liouville boundary bootstrap approach. It is 
remarkable that the above equations are sufficient to reproduce the correlation functions. These 
are much simpler than the standard loop equations, which can be obtained as Ward identities 
in the corresponding matrix model, and which have contact terms associated with degenerate 
world sheets. 

In Liouville theory the difference equations follow from the fusion rules with the lowest 
degenerate boundary fields which shift the boundary parameter. In the microscopic picture, the 
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difference equations are associated with boundary operators, whose geometrical meaning is to 
create an open line on the world sheet starting at the boundary. These operators are associated 
with the lowest degenerate matter fields. It is therefore natural to expect that there is an 
underlying algebraic structure involving operators which are products of Liouville and matter 
degenerate fields. In the case of vanishing boundary cosmological constant such a structure, 
the 'boundary ground ring', has been discovered in |49J . Its generalization to the case of finite 
boundary parameters seems quite straightforward. The ground ring structure should lead to 
difference equations for the boundary n-point function for any n > 2, which is also the case for 
the difference equations derived here in the microscopic approach. 

In this paper we considered only the gaussian field realization of the matter field (without 
screening operators), whose lattice analog is the non-restricted SOS model. The difference 
equations can be of course derived for the 0(n) matrix model and the matrix models associated 
with ADE Dynkin graphs. 
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A Special functions 

• r^(x) function 

The Double Gamma function introduced by Barnes [IZj is defined by: 



d °° 

l0gr 2 (s|wi,W 2 )= ^7 ^2 ( S + niUJl + n 2^2)~* 



dt 

ni,"2=0 , t _Q 



Definition: T b (x) = r ^ /2 ^ b \ 
Functional relations: 



\Z2Trb bx ~^ 

r " (x + 5)= r(fex) Fb{x) > 

r ' (i+1/6) =-r(^r rt(I) - 

rfe(x) is a meromorphic function of x, whose poles are located at x = —nb—mb~ l ,n, m £ N. 
Integral representation convergent for < Rex 



logr fe (x) 
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Sb(x) function 

Definition: S b (x) = r ^Q- x ) 
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Functional relations: 

Sb(x + b) = 2sin(irbx)Sb(x), 
Sb(x + 1/b) = 2sin(irx /b)Sb(x). 

Sb(x) is a meromorphic function of x, whose poles are located at x = —nb—mb~ l ,n, m G N, 
and whose zeros are located at x = Q + nb + mfr -1 , n, m E N. 
Integral representation convergent in the strip < Rex < Q 

sinh(§ - a?)t (Q - 2x) 
2sinh(|)sinh(^) t 



logSb(x) = / - 
Jo 1 



B Functional relations for the boundary three point function 

The LFT boundary three point function satisfies the pentagonal equation |13j : 
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where F, 



ft P2 
/32lfe2 P 4 Pi 

that express the invertib 
blocks. 



is the fusion coefficients (its explicit expression can be found in |48j ) 
e fusion transformation between the s- and t-channel Liouville conformal 



In the case where the boundary field B^' a2 is replaced by the degenerate boundary field 

B b 2 , then the associativity condition (|B.1|) gets replaced by 
— 2 
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and /?32 can take the two values ft ± |. We make the choice ^32 = ft — |. We now use 
the particular values for the boundary three point functions |12| : 
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This leads to the equation 
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A similar equation is obtained by considering the degenerate boundary field B 
(EH): 



in 
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where we introduced C 2 ' CTl ^ 1+2 
Subtracting these two equations gives 



ct(/?i,<7i,cr 3 ). 



^-,(04,03 — |,(Tl)/34 „((T4,o-3 + |,(Tl)/34 

/33-|,/3i ~ ft-fA 



A -I 

04 A 



(cl(/3i, 01,03) - c_(/?i,cri,cr 3 ))F_ + 
From the results of J2] we derive 
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vr 7 (-6 2 ) / T(2 + 6 2 - 2bPx 
and it is well known [2] that the fusion coefficient F |_ is equal to 
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We finally get the following functional relation for the boundary three point function : 
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Another insertion of the spin —6/2 gives: 
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Now we use the relation between the boundary three point function and the fusion matrix 
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where g(Pi, o~ 2 , <J\) is the normalization of the boundary operators computed in 13 ; this 
leads to the following functional relation for the boundary three point function, that pro- 
duces shifts on the boundary conditions only: 

- simrbio-i + a 2 - ft - 6/2) sin7r6(<T 3 + a 2 - 2 - ^(-^A^A^fe 
+ sin^( CTl -a 2 -(5 l + 6/2) sinvr6(a 3 - a 2 - ^c^-W^+b^Ps 
= sin 7r6(<7i + a 3 - fa - 6/2) sin 7t6(2cj 2 - Q)C^ 2 ' ,Tl+f,/2)/33 . (B.9) 



References 

[1] A. Polyakov, "Quantum Geometry Of Bosonic Strings", Phys. Lett. B103 (1981) 207. 

[2] V.Kazakov, "Bilocal regularization of models of random surfaces", Phys. Lett. B150 (1985) 
282; F.David, "Planar diagrams, two-dimensional lattice gravity and surface models", Nucl. 
Phys. B257 (1985) 45; D. Boulatov, V.Kazakov, I.Kostov and A.Migdal, "Analytical and 
numerical study of the model of dynamically triangulated random surfaces", Phys. Lett. 
B157 (1985) 295. 

[3] P. Di Francesco, P. Ginsparg, J. Zinn-Justin, "2D Gravity and Random Matrices", 
Phys.Rept. 254 (1995) 1-133, e-Print Archive: hep-th/9306153. 

[4] P. Ginsparg and G. Moore, "Lectures on 2D gravity and 2D string theory (TASI 1992)", 
|hep-th/93046TT| 

[5] V. Knizhnik, A. Polyakov, A. Zamolodchikov, Mod. Phys. Lett. A3 (1988) 819. 

[6] F. David, Mod. Phys. Lett. A3 1651 (1988); J. Distler, H. Kawai, Nucl. Phys. B321(1988) 
509 . 

[7] M. Goulian and B. Li, Phys. Rev. Lett. 66 (1991), 2051. 
[8] V. Dotsenko, Mod. Phys. Lett. A6(1991), 3601. 

[9] A. Belavin, A. Polyakov, A. Zamolodchikov, "Infinite conformal symmetry in two- 
dimensional quantum field theory", Nucl. Phys. B241, 333-380, (1984). 

[10] H. Dorn, H.J. Otto: Two and three point functions in Liouville theory, Nucl. Phys. B429 
(1994) 375-388 

[11] A.B. Zamolodchikov, ALB. Zamolodchikov: Structure constants and conformal bootstrap 
in Liouville field theory, Nucl. Phys. B477 (1996) 577-605 

[12] V. Fateev, A. Zamolodchikov, Al. Zamolodchikov, Boundary Liouville Field Theory I. 
Boundary State and Boundary Two-point Function, e-Print Archive: hep-th/0001012. 

[13] B. Ponsot, J. Teschner, "Boundary Liouville Field Theory: Boundary three point function", 
Nucl. Phys. B622 (2002) 309, e-Print Archive: hep-th/0110244 

[14] K. Hosomichi, "Bulk-Boundary Propagator in Liouville Theory on a Disc", JHEP 0111 044 
(2001), e-Print Archive: hep-th/0108093, 



50 



[15] I.K. Rostov, "O(n) vector model on a planar random lattice: spectrum of anomalous di- 
mensions", Mod.Phys.Lett.A4 (1989) 217. 

[16] I. Rostov, "The ADE face models on a fluctuating planar lattice", Nucl. Phys. B326 (1989) 
583. 

[17] I.R. Rostov, "Loop amplitudes for nonrational string theories, Phys. Lett. B266 (1991) 317. 

[18] I. Rostov, "St rings with discrete target space", Nucl. Phys. B376 (1992) 539, 
|hep-th/9112Q59 

[19] B. Duplantier, "Conformal Fractal Geometry and Boundary Quantum Gravity", 
|math-p h/0303034l 

[20] B. Nienhuis, "Coulomb gas formulation of 2-d phase transitions", in Phase transitions and 
critical phenomena, Vol. 11, ed. C.C. Domb and J.L. Lebowitz (Academic Press, New York, 
1987), ch.l. 

[21] B. Duplantier and H. Saleur, Phys. Rev. Lett. 63 (1989) 2536; B. Duplantier, Phys. Rev. 
Lett. 64 (1990)493. 

[22] R. Friedrich, W. Werner, "Conformal fields, restriction properties, degenerate representa- 
tions and SLE" , |math.PR/020938"2| M. Bauer, D. Bernard, "Conformal Field Theories of 
Stochastic Loewner Evolutions", e-Print Archive: hep-th/0210015. 

[23] B. Duplantier and I. Rostov, "Conformal spectra of polymers on a random surface", Phys. 
Rev. Lett.61 (1433) 1988. 

[24] V. Razakov, I. Rostov, "Loop Gas Model for Open Strings", Nucl. Phys. B386 (1992) 520. 

[25] I. R. Rostov, "Boundary Correlators in 2D Quantum Gravity: Liouville versus Discrete 
Approach ", Nucl. Phys. B658 (2003) 397, e-Print Archive: hep-th/0212194. 

[26] A. B. Zamolodchikov and A. B. Zamolodchikov, "Liouville field theory on a pseudosphere, 
e-Print Archive: hep-th/0101152. 

[27] J. McGreevy and H. Verlinde, "Strings from tachyons: The c = 1 matrix reloaded", e- 
Print Archive: hep-th/0304224; E. J. Martinec, "The annular report on non-critical string 
theory", e-Print Archive: hep-th/0305148; I. R. Rlebanov, J. Maldacena and N. Seiberg, 
"D-brane decay in two-dimensional string theory", e-Print Archive: hep-th/0305159; J. Mc- 
Greevy, J. Teschner and H. Verlinde, "Classical and quantum D-branes in 2D string the- 
ory", e-Print Archive: hep-th/0305194; S. Alexandrov, V. Razakov and D. Rutasov, "Non- 
Perturbative Effects in Matrix Models and D-branes", e-Print Archive: hep-th/0306177. 

[28] A. M. Polyakov, "Gauge Fields and Space-Time", e-Print Archive: hep-th/0110196 

[29] J. Cardy, "Conformal invariance and surface critical behavior", Nucl. Phys. B240 (1984) 
514; "Boundary conditions, fusion rules and the Verlinde formula", Nucl. Phys. B324 (1989) 
581. 

[30] A. Hashimoto, "Dynamics of Dirichlet-Neumann Open Strings on D-branes" , Nucl. Phys. 
B496 (1997) 243, e-Print Archive: hep-th/9608127. 



51 



[31] I. Affleck, M. Oshikawa and H. Saleur, "Quantum Brownian Motion on a Triangular Lattice 
and c = 2 Boundary Conformal Field Theory" , cond-mat / 0009084 

[32] N. Seiberg, "Notes on Liouville theory and quantum gravity" , Lectures given at the 1990 
Yukawa International Seminar, Common Trends in Mathematics and Quantum Field The- 
ories, May 17-19, 1990 and at the Cargese meeting, Random Surfaces, Quantum Gravity, 
and Strings, May 27 - Jun 2, 1990. Published in Prog.Theor.Phys.Suppl. 102 (1990) 319. 

[33] B.L. Feigin, D.B. Fuchs, Representation of the Virasoro algebra, in: A. M. Vershik, D. P. Zh- 
elobenko (Eds), Representations of Lie groups and related topics, Gordon and Breach, 
London, 1990. 

[34] V.S. Dotsenko and V.A. Fateev, "Four point correlations functions and the operator algebra 
in the two dimensional conformal invariant theories with the central charge c < 1" Nucl. 
Phys. B251 (1985) 691. 

[35] J. Teschner, "On the Liouville three point function", Phys. Lett. B363 (1995) 6. 

[36] J. Teschner, unpublished. 

[37] Al.B. Zamolodchikov, conference on Liouville field theory, Montpellier January 1998, un- 
published. 

[38] I. Runkel, " Boundary structure constants for the A-series Virasoro minimal models" , Nucl. 
Phys. B549 (1999) 563, e-Print Archive: hep-th/9811178. 

[39] R. Behrend, P. A. Pearce, V. Petkova and J.-B. Zuber, "Boundary conditions in rational 
conformal field theories", Nucl. Phys. B579 (2000) 707. 

[40] R. Baxter, Exactly Solved Models in Statistical Mechanics, Academic Press, 1982. 

[41] G. Andrews, R. Baxter and P. Forrester, J. Stat. Phys 35 (1984) 35. 

[42] V. Pasquier, "Two-dimensional critical systems labelled by Dynkin diagrams," Nucl. Phys. 
B285 (1987) 162; J. Phys. 35 (1987) 5707 

[43] V.A.Kazakov, "Ising model on a dynamical planar random lattice: exact solution", Phys. 
Lett. A 119 (1986) 140; D.V.Boulatov and V. Kazakov, "The Ising model on random planar 
lattice: the structure of phase transition and the exact critical exponents" , Phys. Lett. B186 
(1987) 379. 

[44] I.K.Kostov, Gauge-invariant matrix model for the A-D-E closed strings, Phys. Lett. B297 
(1992) 74. 

[45] L.Kadanoff, J. Phys. A 11 (1978)1399; H. Knops, Ann. Phys. 128 (1980) 448; M. den Nijs, 
Phys. Rev. B27 (1983)1674 

[46] B. Duplantier, Phys. Rev. Lett. 57 (1986) 941. 

[47] E. W. Barnes, "The theory of the double gamma function", Phil. Trans. Roy. Soc. A 196 
(1901) 265. 

[48] B. Ponsot, J. Teschner, "Liouville bootstrap via harmonic analysis on a noncompact quan- 
tum group", e-Print Archive: hep-th/9911110 



52 



[49] M. Bershadsky and D. Kutasov, "Scattering of open and closed strings in (1+1)- 
dimensions", Nucl. Phys. B382 (1992) 213, e-Print Archive: hep-th/9204049. 



53 



